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Dynamical and thermal characterizations of excited nuclear systems produced during the colli- 
sions between two heavy ions at intermediate incident energies are presented by means of a review 
of experimental and theoretical work performed in the last two decades. Intensity interferometry, 
applied to both charged particles (light particles and intermediate mass fragments) and to uncharged 
radiation (gamma rays and neutrons) has provided relevant information about the space-time prop- 
erties of nuclear reactions. The volume, lifetime, density and relative chronology of particle emission 
from decaying nuclear sources has been extensively explored and has provided valuable information 
about the dynamics of heavy-ion collisions. Similar correlation techniques applied to coincidences 
between light particles and complex fragments are also presented as a tool to determine the internal 
excitation energy of excited primary fragments as it appears in secondary-decay phenomena. 

PACS numbers: Valid PACS appear here 
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I. INTRODUCTION 

Heavy-ion collisions are the only terrestrial means to 
explore the properties of nuclear matter under extreme 
conditions. In order to extract such nuclear matter prop- 
erties, a clear understanding of^he com|Dkx_^dvnaniics of 
heavy-ion collisions is required 1, 2,*^, 4,'?,'^, 
The detected particles are indeed produced by different 
emission mechanisms and at different stages whose ex- 
perimental identification is challenging. Researchers have 
therefore intensively focused on obtaining a well defined 
dynamical and thermal characterization of particle and 
fragment emitting sources. 

Where and when are fragments produced? What are 
their thermal properties, i.e. excitation energy, internal 
temperature, or spin? At what density do nuclear mul- 
tifragmentation phenomena occur? Can we learn some- 
thing about their link to a liquid-gas phase transition in 
nuclear matter [H IH El El El ? 

In this chapter we will present a review of those re- 
search activities that have been devoted to finding an- 
swers to these questions. We will first focus on the 
space-time characterization of particle emitting sources, 
namely the estimation of their sizes, shapes, densities, 
lifetimes and emission chronology. This task has been 
extensively addressed with intensity-interferometry stud- 
ies by exploring light particle-light particle and IMF- 
IMF (Intermediate-Mass Fragment) correlation func- 
tions. The last section will be devoted to the thermal 
characterization of emitting sources by means of light 
particle-IMF correlation function techniques, providing 
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information about fragment internal excitation energies 
and the relative proportion of the thermal component. 
We will finally conclude with some remarks and perspec- 
tives for future research in this field. 



II. INTENSITY INTERFEROMETRY AND 
LIGHT PARTICLE EMISSION 

The space-time properties of heavy-ion collisions can 
be accessed by intensity interferometry 0,0, 18]. This 
technique was originally introduced in astronomy by 
Hanbury-Brown and Twiss to measure astronomical dis- 
tances, such as the radii of stars and galaxies [l^ [20j . 
The technique was later extended to subatomic physics 
by Goldhaber et al. who studied distributions of K 
mesons in proton- antiproton annihilation processes [2ll| . 
Due to their bosonic nature, two-pion correlation func- 
tions show an enhancement at zero-relative momentum. 
The width of this enhancement provided information 
about the volume of the region emitting pions in the stud- 
ied processes. Pion-pion interferometry plays still today 
a key role in the study of heavy-ion collisions at ultra- 
relativistic energies where it is an important observable 
to investigate the production of the Quark Gluon Plasma 

The use of intensity interferometry in heavy-ion col- 
lisions at intermediate energies is generally character- 
ized by a more complicated scenario, as compared to the 
case of astronomical applications. During a nuclear re- 
action not only photons but several particle species can 
be emitted: neutrons, protons, complex particles or frag- 
ments. These emitted radiations can be either bosons or 
fermions, interacting with one another by means of re- 
pulsive Coulomb and attractive nuclear forces. Another 
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FIG. 1: Two-proton correlation function measured in 
i4n+i97^u collisions at E/A=75 MeV (from Ref. fgg, 30J). 



important complication inherent to heavy-ion collisions 
is represented by the fact that the produced nuclear sys- 
tems live for a very short time ranging between 10^^^ 
and 10^^^ seconds. This situation is very different from 
the case of astronomical objects, where the geometry of 
a static object is studied. Moreover, in nuclear reactions 
different particles can be produced at different times and 
by different sources. Therefore, only a full space-time 
characterization of all these multiple emitting sources can 
improve our understanding of heavy-ion collision dynam- 



A. Measuring two-particle correlation functions 

Given two particles with momenta pi and p2, total 
momentum P = pi + P2 and momentum of relative mo- 
tion q ~ (pi/rni — P2hTi2), the two-particle correlation 

function, 1 + i? (^q, is defined experimentally by the 
following equation: 
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Yi2 (pi,P2) = Ci2 ■ 



1 + R 



Yiipi)-Y2ip2) 

(1) 

fn this equation, Y12 {pi,P2) is the two-particle coinci- 
dence yield while and Yi (pi) and Y2 (^2) are the single- 
particle yields. The normalization constant C12 is com- 
monly determined by the requirement R{q) = a,t large 
relative momentum values, q. fn order to obtain suf- 
ficient statistics, the sums in Eq. (I) are performed 
over all detector and particle energy combinations sat- 
isfying a specific gating condition. Experimental stud- 
ies have therefore focused on two types of observables: 
directionally- gated and angle-averaged correlation func- 
tions. Directionally-gated correlation functions are con- 
structed by selecting particle pairs with specific condi- 
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FIG. 2: Deuteron-alpha correlation function measured in 
Sn-I-Sn collisions at E/A=50 MeV with the LASSA detector 
array |3ll |. 



tions on the relative direction between the relative mo- 
mentum, q, and the total momentum, P 23, 24] . For 
instance, correlation functions with the vector q either 
parallel or perpendicular to the vector P have been ex- 
tensively constructed [2^ . Studies with such direc- 
tional gates generally require high statistics. Directional 
effects might indeed be small and difficult to analyze . 
Alternatively, one can study angle-averaged correlation 
functions by integrating over the relative angle between 
the vectors q and P. The resulting correlation function 
depends only on the magnitude of the relative momen- 
tum, q: 

Y12 {Pi , P2) = Ci2 ■[l + R{q)]-Y,Yi (Pi) • ^2 (A) (2) 

Experimentally, the product of the single yields, Yi (pi) • 
Y2 {P2) in Eqs. (I) and (2), has often been approximated 
with the uncorrelated two-particle yields, Yj^""^" (pi,P2), 
constructed via the so-called event-mixing technique: 
particles I and 2 are taken from two different events 
and the correlation function is calculated as I + R{q) = 
C12 • Yi2/Fi7^° . 

Figures f and 2 show examples of angle-averaged 
proton-proton and deuteron-alpha correlation functions, 
respectively, represented as a function of the relative mo- 
mentum q between the particle pairs |^ [s^, 113 ■ If the 
two emitted particles were totally uncorrelated, the prob- 
ability of detecting them in coincidence would be equal 
to the product of the probabilities of detecting the sin- 
gles, i.e. Y12 {pi,P2) « Yi ipi) ■ Y2 {P2), resulting in a flat 
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correlation function, R{q) = 0, at all g-values. Experi- 
mentally it is easily observed that this is not the case. 
As one detects pairs at small relative momentum, strong 
deviations from unity are observed. These deviations are 
due to quantum statistics and to the so-called final state 
interactions (FSI) 23] . In the case of identical fermions 
(bosons), the relative wave function must respect anti- 
symmetrization (symmetrization) rules that induce mea- 
surable effects in the correlation function at small relative 
momenta |0, |23, |^ . Furthermore, the coincident 
particles can interact with their mutual Coulomb and 
nuclear interaction. The Coulomb repulsion is responsi- 
ble for the anti-correlation at small q- values. The nuclear 
attractive force is responsible for the observed prominent 
peaks both in p — p and d — a correlation functions. 

At large relative momenta, g « oo, the correlation 
hmctions shown in Fig. 1 and 2 appear as flat, R{q) ~ 0, 
indicating the absence of correlations between the coinci- 
dent particles. However, the presence of collective motion 
can generate correlation effects even at large relative mo- 
menta where the correlation function may significantly 
deviate from the limit, R{q) « 0. If collective motion 
exists, the uncorrelated relative momentum distribution, 
constructed by mixing particles from different events, can 
contain additional collective components that do not ex- 
ist in the coincidence spectrum. These additional contri- 
butions may affect the correlation function constructed 
from the ratio of the coincidence and the uncorrelated 
spectra |27ll2l|. 

In general, the shape of the correlation function is sen- 
sitive to the space-time properties of particle emitting 
sources produced during the reaction |23|. In order to 
properly access information about these emitting sources, 
the use of detector arrays with a high angular and energy 
resolution is required. Especially the angular resolution 
plays an important role in determining the exact location 
and the shape of the resonance peaks and in accessing the 
correlation function at very low relative momentum. In 
this respect, position sensitive detectors and silicon strips 
have been quite successful thanks to their capability of 
providing relative angle measurements as small as 0.1°- 
0.3°. 



III. SPACE-TIME PROPERTIES FROM 
TWO-PROTON CORRELATION FUNCTIONS 

Intensity interferometry has extensively been used with 
protons, these particles being abundantly produced at all 
incident energies and easily detected with high resolution. 
Theoretically, the proton-proton correlation function is 
calculated by the so-called Koonin-Pratt equation (KP 
equation) j2J|: 

1 + R{q) = l + J drS{f)-K{f,q) (3) 

The goal of intensity interferometry consists of solving 
Eq. (3): from the measured correlation function on the 
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FIG. 3: Two-proton correlation functions calculated by means 
of Eq. (3) assuming a Gaussian spherically symmetric source 
function, S (r) oc exp {—r^/ro), with source size ro=2.5 (top 
panel) and 5 fm (bottom panel). The different lines cor- 
respond to calculations performed by considering the anti- 
symmetrization and all the final state interactions (solid line), 
only the Coulomb interaction (dashed line) , and the Coulomb 
interaction and the anti-symmetrization of the wave function 
(dotted line). 



l.h.s of Eq. (3) one needs to extract the unknown source 
function, S (r) . This source function is defined as the 
probability of emitting two particles at relative distance 
r, calculated at the time when the last of the two particles 
is emitted. The so-called kernel function, K {r,q), can be 
calculated as K{f,q) = \'^^{r) \ — 1, where ^n(r) is the 
proton-proton scattering wave function 0, |23, 13, 112 ■ 
The kernel contains all the information about the anti- 
symmetrization of the proton-proton wave function, due 
to their Fermionic nature, and the mutual Coulomb and 
nuclear final state interactions (FSI). 

Figure 3 shows correlation functions calculated with 
Eq. (3) by using a Gaussian-shaped source function, 
S (r) (X exp(— r^/rg), with width parameter values of 
rQ—2.5 fm (small source, top panel) and rQ—5 fm (large 
source, bottom panel) jig . Gaussian sources have ex- 
tensively been used in the literature due to their sim- 
plicity. The dashed lines in Fig. 3 correspond to the 
correlation functions obtained for protons interacting 
only with the mutual Coulomb force. This repulsive 
interaction between protons is responsible for the anti- 
correlation at small relative momentum, q<15 MeV/c. 
The dotted line is obtained by adding the two-fermion 
anti-symmetrization in the two-proton wave function, in- 
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FIG. 4: Effects of finite emission lifetimes on two-proton source functions (left side) and directionally-gated correlation functions 
(see text for details). 



ducing a further anti-correlation in the region q= 15-60 
MeV/c, due to the Pauli principle that prevents protons 
from occupying relative momentum states over an inter- 
val of Apx — h/Ax, where Ax is the spread of the proton 
spatial distribution in a certain direction x. Finally, when 
also the final state mutual nuclear interaction is included, 
the p — p correlation functions displays a prominent peak 
at 20 MeV/ c. This peak is due to the s-partial wave of the 
proton-proton scattering problem and strongly depends 
on the volume of the emitting source. Indeed, when a 
large emitting source is used (bottom panel in Fig. 3), 
the peak height reduces significantly. 



The correlation functions displayed in Fig. 3 have 
been determined assuming that the two protons are si- 
multaneously emitted by a source with zero lifetime. 
However, proton emission is known to occur over finite 
timescales (ranging from few tens of fm/c in the case 
of pre-equilibrium emission to thousands of fm/c in sec- 
ondary decay processes). If the particle pairs are not 
emitted simultaneously, the source function is affected by 
a space-time ambiguity and will appear deformed, as it is 
schematically shown in Fig. 4 (left panel) .23, 25] . If tq 
is the actual geometrical source size, the source function 
appears elongated and with a larger size in the direction 
defined by the total momentum vector P. The elonga- 
tion is approximately given by Vt, with V and r being, 
respectively, the average pair velocity and the emission 
lifetime. Due to this deformation of the source function, 
the Pauli suppression effect, described in Fig. 3 (dotted 
lines), plays a key role in determining the shape of the 
correlation function. As it is schematically shown in the 
right panel of Fig. 4, the transverse correlation function, 
constructed with the relative momentum, q, perpendicu- 
lar to the total momentum, P, will undergo a larger Pauli 
suppression as compared to the correlation function con- 
structed by selecting a relative momentum, q, parallel to 
the total momentum, P. This simple qualitative argu- 
ment shows that directionally-gated correlation functions 
allow one to disentangle the space and time information 
hidden in the source function, providing quantitative es- 
timates of both source size and lifetime 12a, U&. . 
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FIG. 5: Data points: directional two-proton correlation func- 
tions measured in Ar-fSc collisions at E/A=80 MeV (from 
Ref. H^). The solid and dashed lines correspond, respec- 
tively, to longitudinal and transversal calculated correlation 
functions. 



A. Volumes and Lifetimes from 
Directionally-Gated correlation functions 

The sensitivity of directional correlation functions to 
finite lifetimes and emission volumes has stimulated an 
extensive research activity in the field of heavy-ion colli- 
sions. 

Figure 5 shows two-proton correlation functions mea- 
sured in central Ar-|-Sc collisions at E/A=80 MeV and 
for proton pairs with total momenta P=400-600 MeV/c 
psf . The filled and open circles correspond, respec- 
tively to longitudinal and transverse correlation func- 
tions. A stronger Pauli suppression in the transverse 
direction is observed, as compared to the longitudinal 
direction. These results are one of the first experimental 
evidences of the predicted lifetime effects in two-proton 
correlation functions. In the same reaction, protons at 
higher total momenta showed no directional dependence 
|33| . consistent with negligible lifetime effects. Indeed, 
these fast protons are dominated by rapidly decaying pre- 
equilibrium sources. The directional correlations in Fig. 
5 were analyzed by using the Koonin-Pratt equation, Eq. 
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(2) , assuming a source function composed by a Gaussian 
spatial distribution and an exponentially decaying time 
profile, S (r, t) cx exp (— r^/2rQ) • exp {—t/r). The simul- 
taneous best-fit of longitudinal and transverse correla- 
tions provided source radii and lifetimes of ro=4.5-4.8 fm 
and T=f0-30 fm/c, respectively. The obtained geometri- 
cal size, ro=4.5-4.8 fm, is comparable to the size of the 
overlapping region in central Ar-|-Sc collisions. The ex- 
tracted short lifetime was interpreted as consistent with 
proton sources dominated by the early pre-equilibrium 
stage of the reaction. These results were also predicted 
by BUU transport models that are expected to well de- 
scribe the dynamical early stages of the reaction js^] . 

Directional correlation functions have also been used 
to characterize proton evaporation from long-lived nu- 
clear systems at relatively low excitation energies j33 . l35j . 
In the study of quasi-compound projectiles produced in 
i29Xe + 27 Al reactions at E/A=31 MeV [33, a life- 
time of the order of r w 1300 fm/c and a source size 
of ro « 2.2 fm were obtained. However, this source size 
appeared small as compared to the size of the evapo- 
rating compound system. The difficulty in explaining 
the results was attributed by the authors to possible 
model dependences in the source parametrization. Sim- 
ilar difficulties have been encountered in the study of 
32s_^natAg reactions at E/A = 22.3 MeV/u Di- 
rectional two-proton correlation functions were used to 
explore both pre-equilibrium emission (fast sources) and 
proton evaporation (slow source) from the produced com- 
pound nuclei. In both cases no significant difference be- 
tween transverse and longitudinal correlation functions 
were observed. While this observation can be expected 
in the case of pre-equilibrium emission, the absence of 
any directional effects for low energy protons evaporated 
from compound nuclei appeared controversial. 

The space-time properties of heavy-ion collisions at rel- 
ativistic incident energies were explored by the FOPI and 
the ALADiN collaborations "2^, '2?, "281121 El El to 
access information about nuclear densities and to study 
the temporal evolution of excited nuclear systems. Fig- 
ure 6 shows directional two-proton correlation functions 
measured in central ^^Zr(^^Ru)-|-^^Zr(^^Ru) collisions at 
E/A=400 MeV A very small Pauli suppression 

in the transverse correlation function (solid point) as 
compared to the longitudinal correlation function (open 
points) is observed. This suppression is consistent with a 
slightly deformed Gaussian source characterized by a neg- 
ligible lifetime and a geometrical size of about ro ~ 5.4 fm 
[2^ . The obtained negligible lifetime can be attributed to 
the presence of a strong collective motion and the conse- 
quent decrease of emission timescales in participant mat- 
ter. Longer emission lifetimes can be expected in the 
decay of spectator matter produced in peripheral colli- 
sions [llj. The left panels in Fig. 7 show transverse 
(closed symbols) and longitudinal (open symbols) two- 
proton correlation functions measured in the decay of tar- 
get spectators produced in Au-|-Au coUisions at E/A=l 
GeV [33. The data from the top to the bottom are 
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FIG. 6: Longitudinal (open dots and solid line) and transver- 
sal (close dots and dotted line) two-proton correlation func- 
tions measured in ^''Zr(^''Ru)-|-'-'''Zr(^'^Ru) central collisions 
at E/A=400 MeV with the FOPI detector 



obtained by imposing different gates on the T^bound ob- 
servable, defined as the sum of the charges of all frag- 
ments with Z>2 and emitted by projectile spectators [lll| . 
Zfcotirid increases with increasing impact parameter and 
decreasing spectator excitation-energy. These directional 
correlation functions provided lifetimes t < 20 fm/c and 
source radii ro « 8 fm, independently of T^bound- From an 
estimate of the amount of nucleons in the decaying spec- 
tator it was possible to deduce densities ranging between 
p/po ~ 0-15 and p/po ~ 0.4. Even in the case of spectator 
decay, the obtained lifetimes appear to be rather short. 
High proton identification thresholds fHlsl could result 
in correlation functions that are dominated by fast parti- 
cles emerging from the short-lived pre-equilibrium emis- 
sion stages of the reaction. On the other hand, other au- 
thors have investigated the intrinsic small sensitivity of 
directional correlation functions to those long lifetimes 
typical of evaporation processes and secondary decays. 
Can one probe the time structure of proton emissions 
when lifetimes of the order of r >10'^-10^ fm/c exist? 
In ReL [s^ it was observed that, if the emission life- 
time is very long, longitudinal and transverse correlation 
functions become undistinguishable and, as a result, a 
quantitative extraction of the lifetime itself becomes dif- 
ficult. 



B. Space-time source extents from angle-averaged 
two-proton correlations 

Angle-averaged correlation functions are constructed 
experimentally by removing any conditions on the an- 
gle between the relative and the total momentum of the 
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FIG. 7: Left panel: Transverse (closed symbols) and longitudinal (open symbols) two-proton correlation functions measured by 
the ALADiN collaboration in the decay of target spectators formed in Au+Au collisions at E/A=l GeV ,3?]- The data from 
the top to the bottom panel correspond to increasing impact parameter, selected by means of the Ztound observable pH |37|. 
Right panel: ratio between longitudinal and transverse correlation functions. 



emitted proton pairs. From this point of view, these ob- 
servables do not require as high statistics as in the case 
of directionally gated correlation functions. In order to 
extract physics information, one needs to use the angle- 
averaged Koonin-Pratt equation 0, [s^, : 

R (q) = 47r J dr-r^ S{r)-K (r, q) (4) 

where K{r, q) is the angle-averaged kernel. The result- 
ing spherically symmetric source function, S{r), contains 
information about its spatial extent and its finite life- 
times folded together in the value of the relative distance, 
r. Therefore, it is not easy to disentangle the space-time 
ambiguity contained in the emitting source. However, 
even under these limited conditions, angle-averaged cor- 
relation functions have provided space-time information 
about particle emission mechanisms. 

In order to solve the angle-averaged Koonin-Pratt 
equation, Eq. (4), and extract information about the 
emitting source function, S{r), different approaches have 
been proposed. In the following subsections we will 
present the main results obtained with Gaussian source 
approaches and imaging techniques. 

1. Gaussian source approaches 

In Fig. 8 angle-averaged two-proton correlation func- 
tions measured in ^^N-|--'^^^Au reactions at E/A=75 MeV 
are shown for three different gates on the total momen- 
tum of detected proton pairs, P=270-390, 450-780 and 



840-1230 McV/c jHES. The dashed hues correspond 
to best-fits of the experimental data with Eq. (4), where 
the source function is assumed to be characterized by 
a Gaussian shape, Sp{r) cx exp (— r^/2rQ) . Gaussian 
source sizes, ro=5.9, 4.2 and 3.4 were obtained, respec- 
tively, for P=270-390, 450-780 and 840-1230 MeV/c. 

The total momentum dependence of two-proton corre- 
lation functions, shown in Fig. 8, has been extensively 
explored in the literature. Figure 9 shows a collection 
of Gaussian source sizes measured in reactions induced 
by different projectiles (^HCjJf N, ^^O and ''"Ar) imping- 
ing on Au and Ag targets [13, • Source sizes are 
represented as a function of the average velocity of the 
coincident proton pairs, Vp — 1/2 (pi +P2) /Mp, normal- 
ized to the beam velocity, w&eom- The observed decrease 
of the size of the emitting source with increasing pro- 
ton velocities has been interpreted as a consequence of 
the coo ling dynamics of the produced nuclear systems 
[T^ ItR l30j | Energetic protons are associated with small 
sources as one would expect in the case of emissions at 
the early stages of the reaction. For less energetic pro- 
tons, however, the extracted source radii increase as one 
would expect at the later stage when the system is cooled 
and expanded. The systematics of two-proton radii with 
the size of the interacting projectile and target nuclei is 
more involved. The data points corresponding to reac- 
tions induced by ^"^N and ^^O projectiles overlap and are 
interpolated by a dashed line. The other two dashed lines 
are obtained by multiplying and dividing the 

1/3 

line by a factor proportional to {Ap / 14) ' , where Ap is 
the mass number of the projectile nuclei. It is clear that 
for high energy protons, Vp/vbeam > 0.5, where fast pre- 
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FIG. 8: Data points: Two-proton angle- averaged correlation 
functions measured in "N-^^^'^Au at E/A=75 MeV [2fll.l3n|. 
The different symbols correspond to different gates in the total 
momentum of proton pairs, as indicated on the figure. The 
dashed lines are calculated assuming Gaussian-shaped source 
functions and correspond to best fits of the height of the peaks 
at 20 MeV/c. 



Long T Pre-equlibrium 



a 4 



o 



Ap=40l 



1^ 



Ap=3 



14N,160 

3He 



0.0 



0.5 



1.0 

V /V 



1.5 



beam 



FIG. 9: Gaussian source sizes obtained from the study of 
reactions induced by ^^N and ^^O projectiles (open symbols), 
^He projectiles (closed circles) and ^"Ar projectiles (closed 
diamonds) on ^^^Au targets (adapted from Refs. [T^ll^'l. 



equilibrium emission plays an important role, the two- 
proton source radius scales approximately with the radius 
of the projectile. This scaling with the radius of projectile 
nuclei is not observed in the case of low energy protons 
that are known to be predominantly emitted by evapo- 
ration and secondary-decay processes. The presence of 



such long-lifetime emissions makes two-proton correla- 
tion functions more difficult to be interpreted. This is 
confirmed by the difficulty in reproducing the shape of 
the experimental data for pairs at low total momenta 
(see Fig. 8). 

The Gaussian radii vq shown in Fig. 9 must be consid- 
ered as upper bounds for the actual source sizes. Indeed, 
since angle-averaged source function fold space and time 
extents together, a small size could correspond both to a 
small emitting volume or to short proton emission times. 
On the other hands, the use of gates on total momen- 
tum or, alternatively, on proton velocities is expected 
to reduce the space-time ambiguities contained in the 
angle-averaged correlation function. Indeed, these gates 
should help to isolate, to some extent, more localized 
emitting sources. Stimulated by these ideas, the CHIC 
collaboration has used statistical event generators to sim- 
ulate total-momentum gated correlation functions and 
estimate emission lifetime and spatial extents in Ni-f-Al 
reactions at E/A=45 MeV j^. The authors explored 
also the existence of a possible overlap of multiple sources 
with different lifetimes (corresponding to pre-equilibrium 
and evaporative emissions) . Due to the selectivity of the 
total momentum gate, it was still possible to estimate 
both a spatial Gaussian source size, i?G=2.7±0.3, and 
an emission lifetime, Tp R:i400±200 fm/c. 

Systematic studies of source sizes in central Ca-|-Ca 
and Au-f Au collisions with varying incident energies have 
been recently performed by the FOPI collaboration [36l| . 
The size of the two-proton source is observed to decrease 
as the incident energy is increased from E/A= 400 MeV 
to E/A = 1500 MeV. In the case of Au-|-Au reactions, 
the source size decreases from 5.0 fm to 4.1 fm. A smaller 
source size change is observed in the case of Ca-|-Ca col- 
lisions. This result is consistent with participant mat- 
ter which is compressed to higher and higher densities 
as one increases the incident energy. This indication is 
quite attractive and further work on the details of these 
two-proton correlation functions can provide useful tools 
to investigate the nuclear equation of state. 



2. Imaging analyses 

The importance of studying the detailed shape of two- 
proton correlation functions was recently addressed by 
introducing an imaging approach to intensity interferom- 
etry Qi l45j . This imaging technique consists of 

extracting the source profile, S{r), by a numerical inver- 
sion of the Koonin-Pratt equation, Eq. (6). No a- priori 
assumptions about the source shape are made, thus con- 
siderably reducing model dependences. Figure 10 shows 
an application of the imaging technique to two-proton 
correlation functions measured in ^''N-l-^^^Au collisions 
at E/A=75 MeV [sl El. The thick hues on the left 
panel of the figure represent the imaged measured cor- 
relation functions. The imaging technique reproduces in 
details the entire shape of the correlation functions. On 
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FIG. 10: Imaging analysis ol data points already shown in Fig. 8 and corresponding to two-proton correlation functions 
measured in ^'^N+^^'^Au reactions at E/A=75 MeV [3C| . The left panel shows imaged correlation functions by thick bands. 
The right panel displays the extracted source function profiles ^3 • 



the right panel of Fig. 10 the extracted imaging pro- 
files are also shown. Contrary to the case of Gaussian 
source analyses shown in Fig. 8, the source size does not 
decrease with increasing proton pair total momentum, 
remaining constant at about 3 fm. This result is very 
different from the systematics shown in Fig. 9, where 
the source size was mostly extracted from the height of 
the peak at 20 MeV/c. This constancy of the source size 
with increasing total momentum of the protons shows 
that the height of the peak at 20 MeV/c does not provide 
unambiguous information about the size of the source 
|4^ . These results have been explained within a sim- 
ple scenario where protons are emitted by the overlap 
of a fast source, corresponding to pre-equilibrium emis- 
sions, and a slow source representing the last stages dom- 
inated by secondary decays of excited fragments OS ■ 
In this schematic model, fast emissions provide a frac- 
tion f of the total proton yields, Yfast — f ■ Y, while 
the slow secondary decay sources provide the remainder, 
ysiow = (1 ^ /) ■ In this limit, the profiles extracted 
from the imaging analysis of the experimental R{q) at 
q>10 MeV/c and shown on the right panel of Fig. 10 are 
strongly dominated by the dynamical source. The size 
of this source is determined by the width or, even bet- 
ter, by the shape of the peak at 20 MeV/c. The detailed 
study of the shape provides also a measure of the rela- 
tive contributions of the fast source, /, and of the slow 
source, 1 — / . This fast/slow contributions represent new 
physics information contained in proton-proton correla- 
tion functions, opening the opportunity to constrain the 
importance of secondary decays. The detailed shape of 



the long-lived portion of the emitting source contributes 
to the source function S{r) at large relative distance val- 
ues, r >10 fm, where the kernel K{q,r) is dominated by 
the Coulomb interaction. These large distances influence 
the correlation function only at very small relative mo- 
menta, q <10 MeV/c, where measurements are difficult 
because of the finite angular resolution typical of the used 
experimental setups [42|. Therefore, the details of these 
slow emitting sources r >10 fm are hard to access. In 
Fig. 10, the observed quasi-Gaussian shapes refer only to 
the fast source. This is demonstrated by the fact that 
their integral over relative distances r is less than unity. 
If the slowly emitting sources could be experimentally ac- 
cessed by measuring very low relative momentum particle 
pairs, it would be possible to image even the tails of the 
source profiles at large relative distances. In that case, 
one would clearly observe how the emitting source can 
deviate significantly from the simplified Gaussian shapes 
assumed in most works on two-proton correlation func- 
tions. Such detailed studies on the shape of two-proton 
correlation functions require high resolution experiments 
that represent a challenge for the future. 



With the introduction of imaging analyses, the sys- 
tematics of proton-proton source sizes shown in Fig. 9 
is expected to change. On the other hand, such high or- 
der shape analyses require high resolution measurements 
that certainly represent a challenging perspective for the 
future. 
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FIG. 11: Two-neutron correlation functions measured in 
Ni+Ni, Ni+Au and Ni+Al collisions at E/A=45 MeV by the 
CHIC collaboration El. 



IV. SPACE TIME CHARACTERIZATION 
FROM UNCHARGED RADIATIONS 

Considerable experimental efforts have also been de- 
voted to the measurement of correlation functions involv- 
ing non-charged radiations (neutrons and photons). 

Figure 11 shows two-neutron correlation functions 
measured in Ni-|-Al, Ni, Au reactions at E/A=45 MeV 
[ill 0, 01 . These measurements are very complicated 
because of cross-talk problems: the interaction and regis- 
tration of the same neutron in two different detectors, can 
substantially affect neutron coincidence measurements 

MM 

The neutron-neutron correlation function is dominated 
by the Fermionic nature of the neutrons and by their 
mutual nuclear final state interaction. Unlike the case 
of proton-proton correlation functions, no Coulomb anti- 
correlation at small relative momenta exists. The strong 
maximum at zero relative momentum is caused by the 
same resonance that brings about the maximum in the 
p — p correlation function and has been used to inves- 
tigate the time-scale of the emission and to study the 
interplay between dynamical and statistical effects [ilf . 
Similarly to what was observed in the case of protons, the 
height of this maximum increases with increasing two- 
neutron total momentum .41. 46, 47,] , indicating reduced 
time scales for the emission of more energetic particles 
emerging from the pre-equilibrium stage. The simulta- 
neous study of p — p and n — n correlation functions 
in Ni-t-Al reactions at E/A=45 MeV by means of the 
Csorgo-Helgesson statistical model (Ref. refer- 
ences therein) has provided comparable source radii of 
i?(5=2.7±0.3 fm for both protons and neutrons. How- 
ever, slightly larger timescales, with lifetimes of the order 
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FIG. 12: Data points: Two-photon correlation function mea- 
sured in Ar-I-Al reactions at E/A=95 MeV with the MEDEA 
array [5^ . The solid lines show calculations of the differ- 
ent effects contributing to the overall shape of the correlation 
function. 



of T„ ~ 600±200 fm/c, appeared to be associated with 
neutron emission, if compared to protons emitted with 
lifetimes of about Tp « 400±200 fm/c 0. 

The enhancement at small relative momenta observed 
in the n-n correlation function in Fig. 11 is due to the 
short-range nuclear final-state interaction. However, if 
neutrons are emitted from a long-lived decaying system, 
the increased average spatial separation between succes- 
sive nucleons will reduce the significance of the final state 
strong interaction. The correlation function is then ex- 
pected to be dominated by quantum statistics effects. 
An attempt to observe such an effect was successfully 
performed by the authors of Ref. in the study of 

compound nuclei formed in ^*0-|-^^Mg reactions at E=60 
amd 71 MeV. The n-n correlation function was con- 
structed with the evaporated neutrons and a dip at small 
relative energies was observed. This anti-correlation is 
generated by the anti-symmetrization of the n-n wave 
function and represents the first experimental evidence of 
a pure Fermionic HBT effect. The n-n correlation length 
is longer than the nuclear interaction length and final 
interaction effects are negligible. 

Two-photon correlation functions have been measured 
with the TAPS and MEDEA arrays 51, 52] to extract 
space-time information about 7 emitting sources. Figure 
12 shows a two-photon correlation function measured in 
36 Ar + 27 Al reactions at E/A=95 MeV [E3. The pure 
intensity interferometry effect is observed only at small 
relative momentum (q<45 MeV/c). The peak at 140 
MeV/c is due to the neutral pion decay into two photons 
while the bump at intermediate q-values is produced by 
badly detected pion decays [s^]. Because of these com- 
plications, two-photon correlation functions were studied 
with sophisticated analysis techniques yielding Gaussian 
source sizes of the order of rp ~2 fm and emission life- 
times of the order r « 4 fm/c. These results seem to 
indicate that 7-7 correlation functions probe mostly the 
dynamical stage of the reaction. Similar conclusions were 
also presented in Ref. [H^l where anticorrelations be- 
tween energetic 7-rays and protons were studied with the 
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MEDEA array in Ar+V collisions at E/A=44 MeV. En- 
ergetic gamma rays can be considered very good probes 
of the spatial and temporal properties of the geomet- 
ric overlap region developed in the very early stages of 
the collision process when nucleon-nucleon collisions are 
more important. 



V. SPACE-TIME CHARACTERIZATION FROM 
LIGHT COMPLEX PARTICLES 

During heavy-ion collisions at intermediate energies a 
large variety of particles and fragments are produced. 
A complete space-time characterization of the reaction 
therefore requires multiple intensity interferometry stud- 
ies extended to several nuclear species. Correlation func- 
tions between light charged particles other than pro- 
tons have indeed been investigated by different authors 
(Ref. 01 s-nd references therein). 

Figure 13 shows a, d — a correlation function measured 
in the reaction at E/A=94 MeV 0. The 

large minimum at small relative momentum is due to 
the mutual Coulomb repulsion and the peaks at relative 
momentum q > 30 MeV/c are due to the nuclear inter- 
action. In particular, the sharp peak around 42 MeV/c 
corresponds to the first excited state of ^Li at E*= 2.186 
MeV and the broad peak around 84 MeV/c stems mainly 
from the resonance at E*= 4.31 MeV with small contri- 
butions from the resonance at E*= 5.65 MeV. The solid 
and dashed lines in Fig. 13 show d — a correlation func- 
tions calculated from Eq. (4) with the kernel function, 
K{q,r), constructed from the d — a relative scattering 
wave function ^^l- The best fit of the integral of the 
first peak at 42 MeV/c provides Gaussian source radii 
ro=4.5 and 3.5 fm for d — a pairs with total energies 
Ea + Ed = 55-120 MeV and 120-200 MeV, respectively. 
As already observed in the case of two-proton correlation 
functions, the height of the peaks increases with increas- 
ing velocity of selected d — a pairs. A closer look at Fig. 
13 clearly shows that the whole shape of the d — a cor- 
relation function is not reproduced. Eq. (4) does not fit 
simultanously the first peak at 42 MeV/c and the second 
peak at 84 MeV/c. The broad peak at 84 MeV/c is bet- 
ter described by using a source radius larger by about 1 
fm as compared to the radius extracted from the best-fit 
of the sharp peak at 42 MeV/c. These difficulties have 
not been resolved yet and require further research. 

Attempts to perform systematic size measurements by 
using the first sharp peak at q=42 MeV/c can be found in 
the literature 17]. For instance. Fig. 14 shows the com- 
parison between source sizes extracted from p correla- 
tion functions (left side) and source sizes extracted from 
the best- fit of the first peak at 42 MeV/c in c? — a correla- 
tion functions (right side) measured in the same reaction 
systems, i.e. at E/A=94 MeV, ^°At+^^'^Au 

at E/A=60 MeV and "N-h^^^Au at E/A=35 MeV (Ref. 
|l7j and references therein). Both p — p and d — a source 
radii decrease with increasing energy of the detected 
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FIG. 13: Data points: Deuteron-alpha correlation functions 
measured in ^"^O-f ^''''Au reactions at E/A=94 MeV for 
different gates on the total kinetic energy of the particle pairs. 
The solid and dashed lines correspond to Gaussian anaiyses 
of the first peak after taking into account the finite angular 
resolution of the used experimental setup. 



pairs. Furthermore, d — a radii appear to be smaller 
than p — p source radii. Even if this result seems to be 
intriguing and requires further research, the difficulties 
in reproducing the overall shape of the d — a correlation 
function described in Fig. 13 indicate that these compar- 
isons between d — a and p — p source radii must be taken 
very cautiously. 

Examples of deuteron-deuteron correlation functions 
measured in Ar-|-Au at E/A=60 MeV are shown in 
Fig. 15 |5^ |53| for different gates in the total energy 
of the deuteron pairs (see different symbols). The anti- 
correlation at small relative momentum, due to the re- 
pulsive Coulomb interaction is the main visible feature 
of these observables. These correlation functions have 
been analyzed with final state interaction models that 
include both the Coulomb and the nuclear interaction 
m m m m m m El . source radH from d-d cor- 
relation functions have been found to be generally larger 
than both p — p and d — a radii 0, |HS Ho, "6?, 62'| . 

More recently d-d correlation functions from 
projectile-like sources produced in Xe-|-Sn reactions at 
E/A=50 MeV El have been measured with the INDRA 
47r array |65i l66[ E3 • By means of a quantum model 
inclu ding both Coulomb and nuclear final state interac- 
tions l6g| , deuteron emission times from the decay of the 
projectile-like fragment at different impact parameters 
were estimated. As the impact parameter decreases from 
peripheral to more central collisions, the emission times 
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FIG. 15: Deuteron-deuteron correlation functions measured 
in '^''Ar+^^^Au collisions at E/A=60 MeV ,5S, 57J. The differ- 
ent symbols correspond to different gates on the total kinetic 
energy of the detected deuteron pairs. The lines show cal- 
culated correlation functions constructed by considering both 
the Coulomb and the nuclear final-state interactions. 
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were observed to decrease from 200 fm/c up to 25 fm/c. 
This result was interpreted as an increasingly important 
contribution of out-of-equilibrium emission as one moves 
towards more and more central events where more ex- 
cited projectile-Uke fragment are produced. A better an- 
gular resolution to measure the d—d correlation function 
at very low relative momentum ,64j and an improved 
knowledge of the d—d final state interaction 58] could 
improve our understanding of deuteron emission mecha- 
nisms. 

Triton-triton correlation functions have been also stud- 
ied by some authors with models including both Coulomb 
and nuclear final-state interactions or with simplified pre- 
scriptions using only the Coulomb repulsion between the 
coincident tritons [53, ISfl, |63> ll^ • The extracted source 
radii are comparable to those extracted from d — d cor- 
relation functions, even if uncertainties in the knowledge 
of the triton-triton nuclear interaction still exist . 

In this subsection we have not discussed the possibility 
that the correlation function might be distorted by the 
Coulomb field of the residual system. Only the mutual 
two-body final state interactions are indeed often used 
in the literature. This approximation is reasonable only 
in the case of correlations between particles with iden- 
tical charge-to-mass ratios causing them to experience 
a similar acceleration in the Coulomb field of the resid- 
ual system. Figure 16 shows proton-alpha correlation 
functions measured in ^"Ar-l-^^^Au reactions at E/A=60 
MeV |5^ |6^ . The open (solid) points refer to particle 
pairs with the proton being faster (slower) than the al- 
pha, Vp > Va {vp < Va)- The correlation function displays 
two peaks at about 15 and 50 MeV/c. The first peak is 
due to the sequential decay of ^B, i.e. — > p-|-^Be 
p+{a+a) with detection of a proton and an a particle 
only. The second peak at 50 MeV/c is related to the un- 
bound ground state of ^Li, which has a lifetime of about 
130 fm/c. 

While the location of the peak at 15 MeV/c is not 
altered by the gate on the proton and alpha particle ve- 
locities, the second peak is observed at slightly higher rel- 
ative momenta if pairs with Vp > Va are selected. Since 
the charge-to-mass ratio of protons is greater than that 
of the alpha particles, the proton will experience a larger 
acceleration in the Coulomb field of the residual nuclear 
system. As a consequence, the relative velocity of the pair 
will be decreased if Vp < Va and increased if Vp > Va, re- 
sulting in a shift of the peak position. After taking these 
distortions into account, p — a correlation functions have 
been analyzed in i^N-^is^Au at E/A=35 MeV H 
and the extracted source radii appear to be comparable 
to the radii obtained from d — d and t — t correlations 
and larger than those extracted from d — a correlation 
functions. 

The existing systematics ll^ shows that two-particle 
source sizes change when different nuclear species are se- 
lected. The obtained differences need further investiga- 
tions but they all point to the fact that different particles 
originate from different emitting stages and sources pro- 
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FIG. 16; Proton-alpha correlation functions measured in 
*°Ar+"^Au collisions at E/A=60 MeV 111,163. The closed 
and open circles are obtained by choosing particle pairs in 
which the proton is, respectively, faster or slower than the 
alpha particle. 

duced during the reaction. The complex scenario emerg- 
ing from these studies calls for higher resolution detec- 
tors and improved theoretical descriptions of final state 
interaction models for two-body correlations in heavy-ion 
collisions. In the next sections we will discuss the work 
that has been done to try to understand time ordering 
in the emission of light particles. The study of emission 
chronology is indeed very important in order to improve 
our understanding of particle production mechanisms in 
heavy-ion collisions. 



VI. PARTICLE EMISSION TIME AND 
CHRONOLOGY 

Particles originating from one specific source are emit- 
ted during a certain time interval. Normally, the rate 
of emitted particles changes during this time interval, 
leading to a specific time distribution for the source. 
If a two-particle correlation function of identical parti- 
cles could be constructed from particles emitted from a 
single source (i.e. if the particle source could be iso- 
lated), and if the spatial distribution would be known, 
then the shape of the correlation function would yield 
information on the shape of the time distribution. Nor- 
mally, neither the spatial nor the temporal distributions 
are known. In this case it is not straightforward to ex- 
tract the spatial and time distributions from the shape 
of the correlation function. To interpret the experimen- 
tal results, source models are often used. These source 
models contain some pre-assumption on the spatial and 
temporal distributions. The shape of these distributions 
can, to some extent, be varied by varying parameters of 
the model (such as radius and temporal width parame- 



ters). Such parameters then represent average emission 
point and average emission time, though it should be 
remembered that such average values are model depen- 
dent. In experimental data the situation is much more 
complex, since it is never possible to completely isolate 
one source from all the other sources present during the 
reaction. The contribution from several sources leads to 
a total time distribution with a complex shape, where the 
effective average emission time will depend on different 
factors. 

What one is usually referring to with emission chronol- 
ogy is a difference in the average emission times between 
two particle types. Though, it is worth noticing that 
the emission times of the two particles may overlap to a 
large extent, and the difference in their respective average 
emission times can be small compared with the width of 
the emission time distributions. Furthermore, if different 
sources contribute, like in peripheral reactions, they con- 
tribute to a complex time distribution of the emitted par- 
ticles. As an example, we can consider two particle types 
emitted from two sources. Different origins for differ- 
ent average emission times can be recognized: one origin 
could be a simple shift of two similar time distributions. 
Another origin could be that the width of the time dis- 
tributions are different, while their shape is quite similar. 
Finally, it is also possible to think that, if more than one 
source is contributing, the relative weight of the sources 
is different for the two particle types, leading to different 
average emission times. In a real reaction, all of the above 
reasons contribute, with a different weight depending for 
example on applied gatesor selected angular ranges. 

Even though the extraction of the emission times and 
sequences is quite involved, there is a wealth of experi- 
mental information that is available. A great advantage 
is to perform simultaneous measurements of both like 
and unlike particle correlation functions. By applying 
different gates, such as polar angle gates, total momen- 
tum gates, directional gates, or velocity gates for unlike 
particles, certain sources are enhanced relative to others. 
For instance, non-equilibrium emission can be enhanced 
by high and intermediate total momentum gates, while 
particles from evaporation and excited fragments can be 
suppressed by shape analyses disregarding the very low 
relative momentum region. Furthermore, single parti- 
cle information, such as energy spectra at different an- 
gles, should be used in the analysis. A systematic study 
can, therefore, to a large extent disentangle the space- 
time characteristics of the contributing sources, thereby 
putting strong constraints on theoretical models. 



A. Emission chronology from like particle 
correlations 

The emission chronology between two particle types 
can, under certain conditions, be determined from like- 
particle correlation functions. If it is valid to assume that 
both particle types are emitted from the same spatial re- 
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gion, a fit with a calculated correlation function, based 
on some source function, to the experimental correlation 
function, will yield an average emission time for each par- 
ticle type. By comparing these average emission times, 
an emission chronology can be inferred. For a review, see 
Rcf. 70] and references therein. The same method can 
also be used to determine differences in emission times 
for the same particle type emitted from different systems 
but with the same spatial region, for instance, different 
neutron emission times from systems similar in size and 
energy content, but with different isospin content (see 
next section). 

The drawback of this method is that the results are 
sensitive to the assumption of emission from the same 
spatial region. Furthermore, the extracted average emis- 
sion times are model dependent, since the average emis- 
sion times depend on the pre-assumption of the shape of 
the (spatial and) temporal distributions assumed by the 
specific source model. 



B. Emission chronology from unlike particle 
correlations 

Model independent information on the emission 
chronology of two particle types can instead be obtained 
from unlike-particle correlation functions. If there is a 
difference in the average emission times, it is possible, 
by suitable gates, to divide the particle pairs into two 
classes with different average distances when the two par- 
ticles " start to interact" . Since the strength of the final 
state interaction depends on the distance between the 
particles, this will lead to a different strength of the cor- 
relation function for the two classes. By comparing these 
correlation functions, the emission chronology can then 
be inferred without any model assumptions. 

Below we present the details of this technique. An im- 
portant assumption for this method to be valid is that 
the particles are emitted independently. Furthermore, 
certain assumptions (which depends on the used gates) 
must be made on the spatial region from which the par- 
ticles are emitted. 



1. Velocity gated correlation functions 

A technique to probe the emission sequence and time 
delay of ejectiles in nuclear reactions was first suggested 
for charged particle pairs, based on the idea that mu- 
tual Coulomb repulsion would be experienced by pairs of 
charged particles emitted with a short time delay. Com- 
parison of the velocity difference spectra with trajectory 
calculations would thus give a measure of the average 
particle emission sequence 0, It^, Ell • 

The technique was extended to any kind of interact- 
ing, non-identical particle pairs in the theoretical study 
of Ref. [l^l- There it was demonstrated that the sen- 
sitivity of the correlation function to the asymmetry of 



the distribution of the relative space-time coordinates of 
the particle emission points can be used to determine 
the differences in the mean emission times by applying 
energy or velocity gates. This effect has been proposed 
for particle pairs such as pd and np [t^I, pTf IZSj], K^K^ 

Velocity gated correlation functions of non-identical 
particles is a very powerful tool to investigate emission 
sequences in nuclear collisions 77] . The basic idea is that, 
if there is an average time difference in the emission times 
of two particles types, there will also be a difference in 
the average distance for particle pairs selected with the 
condition vi > V2 as compared to the pairs selected with 
the complementary condition vi < V2- This is because 
the particle emitted first will, on average, travel a dif- 
ferent distance in the two complementary classes (due to 
the different velocities) before the second particle is emit- 
ted. In particular, the interaction is enhanced for those 
pairs for which the average distance is smaller. This can 
be easily seen if one compares the correlation function 
Ci (g) , gated on pairs vi > V2, with the correlation func- 
tion C2 (?) , gated on pairs vi < f 2 . If particle 2 is emitted 
earlier (later) than particle 1, than the condition vi > V2 
will sample smaller distances (larger correlations) than 
the complementary condition vi < V2. Therefore the 
ratio C1/C2 will show a peak (dip) in the region of rel- 
ative momentum q where there is a correlation and a 
dip (peak) where there is an anti-correlation. Further- 
more the ratio C1/C2 will approach unity both for q 
(since the velocity difference of the two emitted parti- 
cles is negligible) and q ^ 00 (since modifications of the 
two-particle phase space density arising from final state 
interactions are negligible). The exact location of the 
peak and dip in the ratio depends on the source, and in 
particular on the origin of the difference in the average 
emission times. 



2. Experimental results 

At low energies, the analysis of two-particle correla- 
tions and velocity difference spectra has allowed to find 
results consistent with the statistical compound nucleus 
decay for the light products emitted from the reactions 
leo+iog (E,^b=62.5 MeV), and (64 MeV). The 

method could also be used for the determination of fission 
time scales j78i] . 

At intermediate energies, the time sequence of p and d 
has been deduced for the E/A=50 MeV Xe-|-Sn reaction 
studied by the INDRA collaboration. An average emis- 
sion of deuterons ~ 250 fm/c earlier than protons has 
been explained as the result of averaging over a long time 
sequence between pre-equilibrium and thermal emission 
for protons, whereas deuteron emission, resulting mainly 
from hard nucleon-nucleon collisions, is concentrated at 
a few tens of fm/c |6^ . 

First experimental evidence of the emission chronology 
of neutrons and protons deduced from the np correla- 
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FIG. 17; From the E/A=61 MeV Ar+Al reaction measured at KVI, velocity-gated np correlation functions (Cn, filled circles, 
and Cp, open circles) and their ratio (Cn/Cp) Left column: forward measurement. Right column: backward measurement. The 
lines are results from model calculations. From Ref. |79|. 



tion function, was reported in Ref. for the E/A=45 
MeV ^^Ni+^'^Al reaction measured by the CHIC collab- 
oration at LNS, Catania. The experimental results from 
differently gated correlation functions were in qualitative 
agreement with the Koonin-Pratt formalism [2^ . It 
was claimed that for events, selected for high-parallel- 
velocity and high-total-momentum which enhances pro- 
jectile residue and/or intermediate velocity sources, the 
proton is, on average, emitted earlier than the neutron. 

In Ref. the emission time chronology of neutrons, 
protons, and deuterons, was presented for the E/A=61 
MeV ^^Ar-|-^^A1 reaction measured at KVI, Groningen. 
The experimental results showed that the angular and 
total-momentum dependences of the pp and np correla- 
tion functions support a dissipative binary reaction sce- 
nario, where early dynamical emission is followed by sta- 
tistical evaporation. The reverse kinematics utilized in 
the experiment, and fairly high energy thresholds, en- 
hanced the early dynamical emission component in the 
backward measurement. The analysis of velocity- gated 
correlations of non-identical particle pairs yielded de- 
tailed information about the particle emission time se- 
quence. The results from the np backward measurement 



are shown in Fig. 1, right column (from Ref. (7^J). The 
dip in the C„/Cp ratio indicates that neutrons are, on the 
average, emitted earlier than protons. For the forward 
measurement (Fig. 17, left column), the shape of the cor- 
relation function exhibits a correlation at q <40 MeV/c 
and a small anti-correlation at 40< q <100 MeV/c, and 
the pairs with w„ > Vp contributing to C„ interact more 
strongly. This indicates that protons are, on the average, 
emitted earlier than neutrons, a result in agreement with 
that obtained for the Oiab ~ 45° E/A=45 MeV ^SNi-h^^Al 
reaction [77| . The complete sequence of average emission 
times, T, extracted from the E/A=61 MeV ^'^Ar-t-^'^Al 
reaction was the following: for the dynamical emission 
source, t„ < < Tp] for the projectile residue emis- 
sion, Td < Tp < Tn- The interpretation of these results, 
presented in Ref. [t^ , highlights the importance of the 
contribution from the different emission sources. 

The technique of Ref. [T^I has also been applied to 
pairs of non-identical light charged particles produced in 
central collisions of heavy ions in the A=100 mass region 
at a beam energy of 400 MeV/nucleon, measured with 
the FOPI detector system at GSI H^. The difference be- 
tween longitudinal correlation functions with the relative 
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velocity parallel and anti-parallel to the center-of-mass 
velocity of the pair in the central source has allowed the 
extraction of apparent space-time differences of the emis- 
sion of the charged particles. Comparing the correlations 
with results of a final-state interaction model, delivered 
quantitative estimates of these asymmetries. Time de- 
lays as short as 1 fm/c or - alternatively - source radius 
differences of a few tenth fm were resolved. The strong 
collective expansion of the participant zone introduces 
not only an apparent reduction of the source radius but 
also a modification of the emission times. After correct- 
ing for both effects a complete sequence of the space-time 
emission of p, d, t, ^He, a particles was extracted. 

At even higher energy regimes, the above method has 
been used to tackle the problem of the possible observa- 
tion of strangelets in the frame of the distillation process 
following the creation of a quark gluon plasma. In this 
case, strange and anti-strange particles may not be pro- 
duced at the same time in a baryon rich system under 
low bag constant scenarios. Such a prediction has been 
tested using K^K~ correlations in Ref. [sj. 

The knowledge of the emission times and chronology 
is very important in order to extract information on the 
nuclear interaction from complex nuclear collisions. It is, 
however, not straightforward to obtain emission times 
and chronology from a nuclear reaction, and different 
methods need to be combined. Using velocity gates is 
a relativeley new and promising method to, model inde- 
pendently, extract the emission chronology. Examples of 
first experimental results using velocity gates have been 
presented, together with their interpretations. To ex- 
clude other explanations, the method needs to be com- 
bined with other probes and methods to draw the correct 
conclusions. This is a field where further developments 
can be made. 



VII. ISOSPIN EFFECTS AND PERSPECTIVES 
FOR THE ASY-EOS 

The isospin dependence of the nuclear equation of state 
is probably the most uncertain property of neutron-rich 
matter. This property is essential for the understanding 
of extremely asymmetric nuclei and nuclear matter as it 
may occur in the r-process of nucleosynthesis or in neu- 
tron stars. In order to study the isospin-dependent EOS, 
heavy ion collisions with isotope separated beam and/or 
target nuclei can be utilized. In these collisions, ex- 
cited systems are created with varying degree of proton- 
neutron asymmetry and a noticeable i sosp in dependence 
of the decay mechanism is expected [S^. By selecting 
semi-peripheral collisions, the symmetry term of the EOS 
is expected to be probed at low densities, while central 
collisions are expected to be sensitive to the high-density 
dependence. 



A. Theoretical predictions 

Recently, the two-nucleon correlation function has 
been considered as a probe for the density dependence of 
the nuclear symmetry energy jsj, ISJ, |83 ■ In these the- 
oretical studies with an isospin-dependent IBUU trans- 
port model, it was shown that the density dependence of 
the symmetry term of the EOS affects the temporal and 
spatial structure of reaction dynamics by affecting the 
average emission times of neutrons and protons as well 
as their relative emission sequence. For central collisions, 
a stiff EOS causes high momentum neutrons and protons 
to be emitted almost simultaneously, thereby leading to 
strong correlations. A soft EOS delays proton emission, 
which weakens the np correlation |83l l84| . It was found 
that the symmetry energy effect becomes weaker with 
increasing impact parameter and incident energy. Also, 
the strength of the nucleon-nucleon correlation function 
is reduced in collisions of heavier reaction systems as 
a result of larger nucleon emission source |84| . It was 
further found that the momentum dependence of both 
isoscalar nuclear potential and the symmetry potential 
influences significantly the space-time properties of the 
nucleon emission source. Specifically, the momentum de- 
pendence of the nuclear potential reduces the sensitivity 
of two-nucleon correlation functions to the stiffness of the 
nuclear symmetry energy |85|| . 

In spite of the large uncertainties in the symmetry po- 
tential, and in particular in the momentum dependent 
part, the exploratory studies in [11,13113 are very en- 
couraging for using two-particle correlations to study the 
symmetry potential. It is also important to remember 
that the symmetry interaction docs not only influence 
the dynamical emission of particles from the overlap and 
neck-like regions, but also the formation of the residues. 
Therefore also the particles emitted from the residues 
contain information on the symmetry energy. To make 
improvements in the future on the understanding of the 
symmetry potential, it is necessary to have models that 
consistently can describe both the dynamical emission of 
particles and the formation of residues and their subse- 
quent pre-equilibrium and equilibrium emission of parti- 
cles. By applying the experimental conditions and sys- 
tematically comparing calculated and experimental en- 
ergy spectra and gated correlation functions for like and 
unlike particles, it will be possible to obtain hard con- 
straints on the symmetry potential. 



B. Experimental results 

First experimental results have been obtained on two- 
particle correlation functions from systems similar in size, 
but with different isospin content |86|. Small-angle two- 
particle correlation functions with neutrons and protons 
have been obtained from semi-peripheral E/A=61 MeV 
Ar-f 112'^2'^Sn collisions measured at the AGOR cyclotron 
of KYI. 
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The emission from the different sources was enhanced 
or suppressed by introducing angular cuts (intermediate- 
velocity-source emission is enhanced at forward, and 
target-residue-emission at backward angles) and cuts in 
the total momentum (Ptot) of the particle pair, calcu- 
lated in the relevant emission source frame. Figure 18 
(from Ref. |83|) presents the pp correlation function for 
particle pairs selected within the three different gates. 

1. Particles emitted by the intermediate veloc- 
ity source at prompt dynamical emission 
stage, (e.g. first-chance nucleon-nucleon 
collisions), are enhanced by selecting high- 
Ptot pairs in the intermediate velocity source 
frame. For the sample of pp pairs, the angular 
range 30° < < 42° is used for this gate. 

2. Particles emitted by the intermediate veloc- 
ity source at a later dynamical emission stage 
(e.g. neck-emission), are enhanced by selecting 
intermediate-Ptot pairs in the intermediate veloc- 
ity source frame. Again, the angular range 30° < 
9 < 42° is utilized for pp pairs. 

3. Particles emitted by the target residue are en- 
hanced by selecting low-Pfot pairs in the target 



residue frame. The angular range 54° < < 120° 
is applied to both pp and np pairs. The neutron 
energy threshold is set to 8 MeV for this gate. 

One can notice that the height of the peak at q«20 
MeV/c is progressively reduced going from gate 1 to gate 
3, indicating an increase in the particle emission time. 
Figure 19 presents the corresponding results for the np 
correlation function. By comparing the results for the 
two Sn targets, one can see an isospin effect, particu- 
larly in gate 1 for pp and in gates 2 and 3 for np, where 
the height of the correlation function is larger for the 
more neutron-rich target. This indicates a shorter aver- 
age emission time for this system. 

For the interpretation of the correlation data, it is im- 
portant to note that the correlation function depends on 
the space-time extent of the emitting source. From the 
size of the source, a stronger correlation is expected for 
the smaller '^^Ar-l-^^^Sn system, an effect expected also 
because of the larger excitation energy per particle avail- 
able for this system (yielding a shorter emission time). 
On the other hand, the change in neutron number im- 
plies a different symmetry energy which also affects the 
n (and p) emission times. Neutrons are expected to be 
emitted faster in the neutron-rich system, which would 
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lead to an enhancement of the correlation strength for 
'^^Ar+^^^'Sn. Thus, the net influence on the correlation 
function is not easily predictable, both due to the un- 
certainty in the symmetry energy and to the presence of 
more than one source of emission. The stronger np cor- 
relation observed for the larger Ar-f ^^^Sn system in the 
low-total-momentum gate may indicate that the more 
asymmetric system generates a more asymmetric and ex- 
cited target residue that, consequently, decays on a faster 
time scale. 

More insight into these results has been gained by per- 
forming an analysis of the particle emission time sequence 
in Ref. In all studied angle and total- momentum 

gates, it was found that neutrons are, on average, emit- 
ted earlier than protons. Furthermore, the shorter np 
emission time scale for the Ar-f ^^"'Sn system results from 
a faster emission of the neutrons, as compared to the 
Ar-f-'^^^Sn system. This is particularly true for the par- 
ticles emitted from the target residue, indicating that 
the residues in the two reactions were formed differently 
due to the symmetry interaction. Further experimental 
results of particle emission sequence involving deuterons 
indicate that for the Ar-|-^^^Sn system neutrons are emit- 
ted slightly earlier than deuterons, again a result point- 
ing to a faster neutron emission for the more neutron rich 
system. Deuterons, being formed mainly by coalescence, 
appear to have emission times that fall in between that of 
neutrons and protons. No sizeable isospin effects in the 
emission sequence of deuterons and protons were found 
for the above mentioned systems. 



VIII. ACCESSING THE SPACE-TIME 
PROPERTIES AT FREEZE-OUT 

IMF-IMF correlation functions are expected to provide 
the space-time properties of nuclear systems produced in 
nuclear reactions at the time when they approach the 
freeze-out stage. This can offer a unique opportunity to 
better understand the mechanisms of multifragmentation 
phenomena and their possible links to a liquid-gas phase 
transition in nuclear matter. 

IMF-IMF correlation functions are usually constructed 
by combining together fragments with different charges 
and masses into a single correlation function. This is 
accomplished by sorting all IMF-IMF pairs with respect 
to the so-called reduced velocity, Vred = Vrei / \/ Z-\ + Zy , 
where Z-] and are the charges of the two IMFs [S^, 
H m m m 11 il m . Typical measured correlation 
functions are shown in Figs. 20-24 and 26. 

Theoretically, final state interaction models including 
only the Coulomb repulsion have often been used to an- 
alyze IMF-IMF correlation functions. Indeed, fragment 
spatial separations are on average expected to be larger 
than the range of the nuclear force. This approxima- 
tion needs certainly to be verified within the context of 
more quantitative models. However, the most promi- 
nent feature observed in IMF-IMF correlation functions 



+ 75 



- Koonin-Pratt Eq. 


■ ■ 1 ■ • ■ ■ 1 ' ' ' ' 1 


- 


Rs=12fm, T=15MeV 






Carbon fragments 

. ' 






' / ^ 




- 

■ 


— ^ / i' ■' 


Classical r (fm/c) 


- 


/ F : 


50 






100 




' / 


200 






500 






° Quantum 






, 1 , . . . 1 . . , , 





5 10 15 20 25 30 



FIG. 20: IMF-IMF correlation functions calculated for differ- 
ent emitting source lifetimes ,89.1 . The lines and dots corre- 
spond, respectively, to classical and quantum calculations. 



is represented by the strong Coulomb anti-correlation at 
small relative velocities (see Figs. 20-24,26). Possible 
resonances corresponding to the decay of heavier frag- 
ments either occur at higher relative velocities not easily 
accessed experimentally or require higher energy and an- 
gular resolution in order to be clearly resolved. Such 
limitations limit the study of IMF-IMF correlation func- 
tions by means of simplified Coulomb interaction mod- 
els. The IMF-IMF Kernel function in Eq. (4) has been 
calculated either from the relative Coulomb wave func- 
tion or by using a classical Coulomb treatment leading to 

1/2 

the simple expression, K{r,q) = {1 — r/rc) — 1, with 
Tc = 2/iZiZ2e^/(7^ and Zi and Z2 being the charges of 
the two coincident fragments |89| |. The emitting source 
function is commonly assumed to have a spherical shape 
with radius Rs and it is assumed to be characterized by 
a temperature T. Fragments are then emitted according 
to an exponential time profile, P{t) oc exp(— i/r). The 
radius, Rs, and the lifetime, r, are deduced as free pa- 
rameters from the best fit of the experimental correlation 
functions, thus providing the space-time characterization 
of IMF emission. Figure 20 shows correlation functions 
calculated with sources having emission lifetimes between 
50 and 500 fm/c. It can be observed that the shape 
of the correlation function at small reduced velocities 
{Vred <30-10~'^c) is strongly affected by fragment emis- 
sion times. The calculations in Fig. 20 are performed by 
using both the classical (lines) and the quantum (dots) 
two-body kernel function in Eq. (4), providing virtually 
identical results 89]. The width of the Coulomb anti- 
correlation at small reduced velocities is used to extract 
IMF emission times. 

Alternative 3-body approaches to IMF-IMF correla- 
tion functions are based on trajectory calculations where 
two IMFs are emitted sequentially with some average 
time delay, and their motion is propagated to the de- 
tectors under the influence of both their mutual flnal 



18 




5 10 15 20 25 30 



FIG. 21: Two-body (lines) and three-body (dots) trajectory 
calculations for IMF pairs with total momenta per nucleon, 
P/A<110 MeV/c (top panel) and P/A>110 MeV/c (bottom 
panel). The calculations are performed for different values of 
the IMF emission lifetimes. 



state Coulomb interaction and the repulsion induced by 
the Coulomb field of the residual system [ool loil l92j| . 
The correlation function is then deduced from Eq. (2) 
applied to the simulated coincidence and single-particle 
yields |93,0,0|. Figure 21 is taken from Ref. ^ and 
shows the comparison between two-body classical calcu- 
lations (lines) and three-body Coulomb trajectories cal- 
culations (points) for a source having charge Zs = 93, ra- 
dius Rs = 12 fm, temperature T = 15 MeV and emission 
lifetimes varying between t — 50 fm/c and r = 500 fm/c. 
The upper and lower panels show the results for carbon 
pairs emitted with total momenta, P/A < 110 MeV/c 
and P/A > 110 MeV/c, respectively. The use of a very 
large mass number for the source (As = 10000) ensures 
that recoil effects are not taken into account and only the 
effects of the residual Coulomb field are explored. For low 
momentum gates, fragments are emitted with small ini- 
tial velocities and distortions in the Coulomb field of the 
source are large. Indeed, significant differences between 
three- and two-body calculations exist up to r=200 fm/c. 
For longer emission lifetimes, the two-body approach still 
remains valid. As it can be easily expected, in the case 
of fragment pairs at high total momentum (lower panel) , 
fragments move so quickly that the effects induced by the 
Coulomb field of the residual system acting as a third 
body can be neglected. The two-body approach is found 
to be still reasonable if the emission times are longer 
(r > 200 fm/c). For shorter emission times, the disagree- 
ment is substantial and three-body correlations seem to 
be important, even if the effect is less pronounced than 
in the case of low momentum fragments (upper panel in 



Fig. 21). Three-body trajectories and even more realistic 
N-body Coulomb trajectory calculations have extensively 
been used in the literature [qoIIqIIoIIo^. 

More recently IMF-IMF correlation functions have 
also been compared to predictions of microscopic mod- 
els where both the Coulomb and the nuclear interactions 
are taken into account [o^, HH . Microscopic models pro- 
vide a more realistic approach and represent certainly 
a promising opportunity for the future in order to have 
a clear and unambiguous space-time characterization of 
complex fragment emission mechanisms. 



IX. CHARACTERIZING 
MULTIFRAGMENTATION PHENOMENA 

Understanding whether multifragmentation results 
from a sequence of binary splittings or rather from a si- 
multaneous break-up of an excited nuclear system has 
certainly represented one of the main questions raised in 
the last two decades. A sequential binary splitting would 
correspond to cluster emission from the surface of an ex- 
cited source (similar to fission). This process is associ- 
ated with long emission times of the order of 10^^"-10^^^ 
sec, necessary for shape deformation. In contrast, if mul- 
tifragmentation corresponds to a simultaneous breakup 
of nuclear matter, the system is expected to fall apart 
over shorter times (10~^^-10~^'^ sec), comparable to the 
timescales involved in the growth of density fluctuations 
in the spinodal instability region of the nuclear phase 
diagram. In the following part of this section we will 
present some of the main results that can be found in 
the literature on the study of fragment emission times in 
heavy-ion collisions. 

One of the earliest studies of fragment-fragment corre- 
lation functions was presented in Ref. (9Qi | where the au- 
thors have studied ^^O-l-^^'^Au and ^'^O+'^^'^Ag reactions 
at E/A=84 MeV. The IMF-IMF correlation function in 
events where two or more IMFs were observed (Fig. 22, 
bottom panel), was compared to the correlation func- 
tion IMF-HR (IMF-Heavy Residue, Fig. 22 top panel) 
in events where only one IMF and one heavy fragment 
were observed. The lines in the bottom panel of Fig. 22 
show the results of three-body Coulomb trajectory cal- 
culations. IMF emission times of the order of t=1000 
fm/c were obtained in both classes of events selected by 
the authors. These results suggest that the production 
of more than two IMFs in the studied reactions is char- 
acterized by a mechanism similar to the evaporation of 
one IMF by a heavy fragment. This result supports a 
sequential binary evaporation mechanism as responsible 
for IMF emission. 

A similar analysis was performed also in Ref. 
where incomplete fusion reactions with the production 
of two or three heavy (Z>10) fragments in ^^Ne-|-^^^Au 
collisions at 60 MeV/u were investigated. The analysis 
of relative velocities and angles in the center-of-mass of 
the coincident fragments was found to favour a sequen- 
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FIG. 22: Bottom panel: The data points show the IMF- 
IMF correlation function measured in ^*0-|-""'Ag reactions at 
E/A=84 MeV |92j in events where at least two IMS fragments 
are detected. The lines correspond to calculations performed 
with three-body Coulomb trajectories assuming different IMF 
emission lifetimes. The top panel refers to events where only 
one IMF and one heavy residues (HR) is observed and shows 
their correlation function. 



tial emission of fragments with short time steps between 
two consecutive binary decays [9^ . 

The data points in Fig. 23 show Be-Be, B-B and C-C 
correlation functions measured in central Ar-|-Au colli- 
sions at E/A=35 MeV at polar angles 12° < 9 < 35° 
[gif. This angular region allowed the authors to select 
IMFs produced from incomplete fusion mechanisms (Ref. 
[9ll | and references therein). The curves in Fig. 23 repre- 
sent calculations performed with Eq. (4) using a classi- 
cal approximation for the two-fragment Coulomb kernel 
function 91]. IMFs were assumed to be emitted from 
a spherical source of radius i?s=12 fm and with emis- 
sion times distributed according to an exponential law, 
P{t) oc exp(— t/r). The comparison to experimental 
data suggests that fragments are emitted with time-scales 
of about T « 100-200 fm/c. Comparable emission times 
were obtained with three-body Coulomb trajectory calcu- 
lations where the IMFs move towards the detectors under 
the influence of both their mutual final-state Coulomb in- 
teraction and the repulsion induced by the Coulomb field 
of the residual system 91]. The obtained emission times 
are shorter than the timescales characteristic of evapo- 
ration processes from compound nuclei. An analysis of 
total-momentum-gated IMF-IMF correlation functions in 
the same set of data showed that high energy IMFs are 
produced with even shorter emission times, of the order of 
50 fm/c, thus indicating that fragment emission in central 
collisions begins in early stages of the reaction and con- 
tinues throughout the later equilibrium stages |9]].l93]. 
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FIG. 23: Data points: beryllium-beryllium (top panel), 
boron-boron (middle panel) and carbon-carbon (bottom 
panel) correlation functions measured in Au-|-Au central col- 
lisions at E/A=35 MeV foU. Lines: calculations performed 
with Eq. (4) using a classical approzimation to the ker- 
nel function and assuming different emission lifetimes and a 
spherical source size of 12 fm. 



The reaction system '^^Ar-l-^^^Au has also been stud- 
ied at incident energies, E/A=50 MeV, using an N-body 
Coulomb approach |9J] . Angle-averaged correlation func- 
tions contain space-time ambiguities that are difficult 
to resolve: the wide minimum of the correlation func- 
tion at small reduced velocities cannot be associated 
with a unique combination of source size, i?s, and emis- 
sion time, T. A long-lifetime emission produces emitting 
sources elongated in the longitudinal direction defined by 
the total momentum vector, P. Similarly to the case of 
p ~ p correlation functions already discussed previously, 
IMF-IMF directional correlation functions were used to 
extract emission times, r = 50 fm/c, from the surface of 
a dilute source having a density of about p/po ~ 0.4 [9^ . 

Central collisions between nearly symmetric systems, 
«6Kr-|-93Nb, at E/A=50 MeV were studied in Ref. ^ 
with the goal of disentangling instantaneous from sequen- 
tial multifragmentation break-up scenarios. Figure 24 
shows correlation functions constructed with IMF frag- 
ments having charges 4< Zjmf <9 detected at polar 
angles 5° < 9iab < 25° (left panel) and 25° < 9iab < 50° 
(right panel). N-body Coulomb trajectory calculations 
(lines in Fig. 24) indicated that the data are consistent 
with very short emission times, r <100 fm/c. The au- 
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FIG. 24: Data points: IMF-IMF correlation functions measured in central collisions ^^Kr+^^Nb at E/A=50 MeV at polar 
angles 5° < < 25° (left panel) and 25° < < 50° [9R| . N-Coulomb trajectory calculations are shown by lines and correspond 
to different emission lifetimes for the detected fragments. 



thors of Ref. have also observed that the study of 
kinetic energy gated IMF-IMF correlation functions is 
sensitive to the details of the fragment emission topology 
(surface/volume emission processes) in the initial state. 
However, no definitive quantitative results could be de- 
duced because of the unknown sensitivity to other corre- 
lations and conservation laws in the initial state [95j | . 

The perspective of extracting emission times from 
IMF-IMF correlation functions seems to be quite attrac- 
tive in order to better understand even the evolution 
of emission mechanisms with the excitation energy de- 
posited into the excited system undergoing the multifrag- 
mcnt decay. The study of such evolution with the excita- 
tion energy has been performed both in central collisions 
at intermediate energies and in peripheral collisions at 
relativistic energies. In the next section we will present 
the results obtained from IMF emission-time measure- 
ments in central collisions where the excitation energy 
is controlled by the incident beam energy. In the fol- 
lowing section we will show similar emission-time studies 
extended to the decay of excited target spectators af- 
ter bombardment by relativistic-energy light probes. In 
this case, the decay of the system is mostly governed 
by the deposited excitation energy, without strong con- 
tributions from collective motion that represent a rele- 
vant phenomenon in the case of central collisions. These 
two sections will allow us to have an idea about how the 
emission times of IMFs are correlated with the degree of 
excitation of the decaying system. 



A. Emission timescales in central collisions: 
sensitivity to the incident energy 

Mean emission lifetimes for multifragment final states 
produced in Kr-|-Nb reactions at incident energies 
E/A=35, 45, 55, 65 and 75 MeV/nucleon were studied 
in Ref. ;9?|. The width of the Coulomb dip in the IMF- 
IMF correlation function at small reduced velocities was 



observed to increase as the bombarding energy increases 
from E/A=35 to E/A=55 MeV. This result is consistent 
with a reduction of IMF emission times as the violence of 
the collision is increased. The measured IMF-IMF cor- 
relation functions were compared to classical three-body 
Coulomb t rajec tory calculations performed with the code 
MENEKA 'lOO'l and the extracted emission times are rep- 
resented in Fig. 25 as a function of the incident energy. 
Long emission times, r « 400 fm/c, are observed at lower 
incident energies, E/A=35 MeV. As the incident energy 
is increased, the emission of IMFs is observed to occur 
over shorter timescales, until a saturation at a value of 
T « 100 fm/c is observed at E/A=55 MeV [lO^. This re- 
sult is consistent with an evolution of the fragment emis- 
sion mechanism from long-lived sequential processes, at 
lower incident energies, to more a simultaneous scenario, 
at higher incident energies. The extracted minimum life- 
time, r w 100 fm/c, is comparable to the timescales of 
growing density fluctuations for nuclear matter in the 
low-density instability regions of the nuclear phase dia- 
gram. 

The results shown in Fig. 25 are based on the assump- 
tion that the fragments are emitted from a system with 
a single freeze-out condition well localized in time. This 
scenario might be too simplified. Figure 26 shows IMF- 
IMF correlation functions measured in Kr-|-Au reactions 
at E/A=55 MeV for three different values of Vmin, de- 
fined as the minimum velocity of the less energetic frag- 
ment of each pair. It is clearly seen that faster IMFs 
exhibit a wider Coulomb hole, indicating shorter emis- 
sion time-scales. Si milar results were observed also at 
E/A=35 and 70 MeV ^0^. The dependence of the emis- 
sion lifetime on the velocity of the emitted fragments can 
be considered as an indication that a single freeze-out 
condition for fragment emission might not exist. This 
indication was confirmed by the calculations performed 
with the Expanding Evaporating Source model 102] that 
predicted decrea sing emission times for fragments with 
higher velocities |lOlj |. Models based on a simultaneous 
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FIG. 25: IMF emission lifetimes extracted in the study of cen- 
tral Kr+Nb collisions at incident energies between E/A=35 
and 75 MeV 
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FIG. 26: IMF-IMF correlation functions measured in Kr-f Au 
central collisions at E/A=55 MeV and corresponding to dif- 
ferent gates on the velocity of the slower of the two coincident 
fragments |101^| . 



multifragmentation scenario were not capable of repro- 
ducing such evolutionary emission pattern. A similar 
evolutionary scenario in multifragment emission is also 
supported by the extraction of different emission times 
for different nuclear species ^103j. For instance, in the 
study of Kr-|-Au collisions at E/A— 70 MeV, emission 
times of about 50 fm/c and 100 fm/c were foun d, re- 
spectively, for carbon and beryllium isotopes |l03j . 

The described results obtained in the study of central 
collisions have provided evidence for short emission times 
at higher energies. While this finding suggests a simul- 
taneous scenario for multifragmentation, velocity-gated 
and Z-gated correlation studies suggest that evolution- 
ary decay processes cannot be excluded |l0lL Il02l llO.'^ . 
The definition of a freeze-out stage might be considered 
more involved than expected, requiring further experi- 
mental and theoretical investigations. In this respect, it 
would be important to compare experimental correlation 
functions to predictions of microscopic models that de- 
scribe the whole dynamical evolution of multifragment- 




E*/A (MeV) 

FIG. 27: IMF emission lifetimes as a function of the deposited 
excitation energy measured in target spectators decay induced 
by the bombardment tt" and p beams at incident momentum 
of 8.0 and 10.2 GeV/c !9d|. 



ing nuclear systems. Similar studies have recently been 
performed ^98, 104J by comparing IMF-IMF correlation 
functions for central Xe -I- Sn and Gd -I- U collisions at 
E/A=32 MeV, measured with the INDRA detector, to 
simulations performed with the BOB model (Brownian 
One Body dynamics) based on the BNV (Boltzmann- 
Nordheim-Vlasov) approach [o^ llOSl Il06j . The authors 
have studied higher order correlations by selecting differ- 
ent decay channels based on fragment charges and spe- 
cific event topologies [osL Il04j . The conclusions of these 
studies seems to support freeze-out times of the order of 
200-240 fm/c and fragment spatial distributions of the 
order of 3-4 Vq in Xe-|-Sn and 8 Vq in Gd+U reactions, 
with Vo being the volume of the source at normal density, 
Vo = (4V3)(1.2)3^tot fm3. 

The study of central collisions is generally complicated 
by the presence of collective motion. The explosive na- 
ture of such systems can decrease the effective IMF emis- 
sion times. It is, therefore, not easy to correlate the IMF 
emission times shown in Fig. 25 with the deposited ex- 
citation energy. 

In the next subsection we will present measurements of 
IMF emission times in target spectators where estimates 
of the actual thermal deposited excitation energy were 
made. Such studies are less dependent on the dynamical 
effects induced by collective motion and can provide bet- 
ter links to the occurrence of a liquid-gas phase transition 
in nuclear matter driven by thermal effects. 



B. Evolution of emission times with excitation 
energy in target spectators 

Fragmentation phenomena of target-spectator frag- 
mentation induced by light probes at relativistic ener- 
gies are one of the best tools to investigate thermally 
driven phase transitions. In the study of 7r~,p-|-Au at 
8.0, 8.2, 9.2 and 10.2 GeV/c, IMF-IMF correlation func- 
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tions from the decay of Au target spectators were studied 
|96j . The evolution of IMF emission times with the exci- 
tation energy per nucleon is represented in Fig. 27. As 
the deposited excitation energy increases, emission hfe- 
times decrease from r « 500 fm/c at excitation energies 
E*/A 2.5 MeV to a saturating value of about t « 20-50 
fm/c for excitation energies above 5 MeV/nucleon. These 
results indicate a transition from a surface evaporation- 
like emission at low excitation energies towards a bulk 
simultaneous multifragmentation scenario above excita- 
tion energies of the order of E*/A=5 MeV [o^. Fur- 
thermore, the extracted emission times seem to be com- 
parable with timescales of thermodynamical fluctuations 
leading to liquid-gas phase transitions in nuclear matter. 
The decreasing emission times should be related to the in- 
creasing thermal excitation energy deposited into the sys- 
tem. Contrary to the case of central collisions presented 
in the previous subsection, only a small collective-motion 
component caused by the thermal expansion of the sys- 
tem should exist and the dynamics of the decay is mostly 
dominated by thermodynamical aspects. 



X. INTERNAL EXCITATION ENERGY OF 
COMPLEX FRAGMENTS 

Light-particle-fragment cor- 
relation functions have been 
used to perform a thermal characterization of the 
emitting source by providing experimental information 
on the size and t he in ternal excitation energy of the 
primary fragments [TotI EQE llOQl [UOl llll| . 

Several statistical and dynamical models predict very 
different internal excitation energies of the primary frag- 
ments. Tli^qriantum molecular dynamic (QMD) model 
0, III2I III3I Ill4j and the microcanonical metropolis 
Monte Carlo (MMMC) model |TT5l | predict a rather 
low excitation energy for primary fragments. Antisym- 
metrized molecular dynamics 116|. ITitI UJM (AM D), 
stochastic mean field simulations [qTLIT^ flO ^ | ll9l . the 
statistical multifragmentation model |l20l [1211 (SMM) 
and microcanonical multifragmentation models 122, 123, 
Il24j predict moderately hot primary fragments. More- 
over, the mechanism responsible for the formation of frag- 
ments differs from a model to another. In dynamical 
models, the formation of the fragments and their exci- 
tation energy depend not only on the collision dynamics 
but also on the pr ocedure emplo yed for performing the 
cluster recognition 1251, 121 [123. On the other hand, in 
the case of statistical models, the excitation energy of the 
fragments is an assumption of the model itself used to cal- 
culate the statistical weights of the partitions. Therefore, 
experimental estimates of the importance of secondary 
decays and of the size of the primary fragments can pro- 
vide very significant constraints and crucial tests for the 
theoretical models. A direct measurement of these quan- 
tities would be also interesting for various aspects: the 
extraction of certain physical information, such as the 



rate of statistically to dynamically emitted particles or 
the caloric heat capacity 0, , would be less depen- 
dent on model assumptions. 

Previous studies have shown that it is possible to ex- 
tract the intrinsic properties of the fragments indepen- 
dently of the mechanism of their formation. These stud- 
ies are based on the measurement of relative- velocity cor- 
relation functions between fragments and light charged 
particles (LCP). The aim of this section is to give a re- 
view of the work dedicated to the determination of the 
sizes and excitation energies of primary fragments. An 
excitation function of these quantities will be given for 
Xe -I- Sn s ystein for the i ncide nt energy range E/A = 
30-50 MeV [lol [lo3 [ni |m| . We will present results 
on th e stu dy of central Kr-|-Nb collisions at E/A — 45 
MeV |l08j and of quasi projectiles produc ed in periph- 
eral Xe+Sn colHsions at E/A = 100 MeV ^2^. These 
studies were conducted by using either the Washington 
University Dwarf Ball/ Wall array or the in INDRA 
multidetector [6^ |6^ ■ It is important to note that 
the same data can provide information on the space-time 
extend of the outgoing fragments, but in this section we 
will concentrate only in the measurement of their excita- 
tion energy. 



A. Reconstruction of the primary fragments 

At intermediate energies, there are at least two ex- 
treme mechanisms for the production of light particles, 
(i) in nucleon-nucleon collisions and (ii) by statistical 
evaporation from excited sources. The former occurs dur- 
ing the first stage of the collision when the dynamics play 
an important role. This process is called direct or pre- 
equilibrium emission. The second process is supposed to 
occur at later stages of the reaction as a secondary emis- 
sion from excited primary fragments. Between these two 
stages, continuous emission may occur from other sources 
that are difficult to define. In multifragmentation events, 
it is possible to isolate the secondary statistical compo- 
nent if the fragments formed are not too excited, so that 
the time scale associated with their decay is much larger 
than the time scale of their production. 

In this context the primary fragment excitation en- 
ergy can be deduced from the multiplicity of its associ- 
ated evaporated LCPs. Figure 28, taken from Ref. |l08j . 
shows the projection of the relative velocity between the 
heaviest fragment and the protons detected in coinci- 
dence and in the case of Kr-t-Nb collisions at E/A = 45 
MeV. The components V|| and V^ are the projections of 
the relative velocity onto the axis representing the frag- 
ment direction in the center of mass frame and into a 
plane perpendicular to that axis. The top panel corre- 
sponds to raw coincidence data, the middle panel corre- 
sponds to the background distribution and the bottom 
panel displays the result of a subtraction of this back- 
ground from the total LCP emission. This bottom panel 
clearly shows a ring surrounding the IMF location in ve- 
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FIG. 28: Reaction system: Kr + Nb at E/A = 45 MeV. Paral- 
lel velocity-transverse velocity invariant diagram for protons 
in the center of mass of the heaviest fragment detected in co- 
incidence in each event. Top panel: correlated events; Middle 
panel: uncorrelated events; Bottom pan el: difference between 
correlated and uncorrelated events |10^ . 
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FIG. 29: (a) Phosphorus-alpha correlation measured in cen- 
tral collisions of Xe -I- Sn at E/A = 32 MeV. (b) Difference 
function, (c) Velocity spectrum of a particles in the center of 
mass of the phosphorus fragments, obtained by subtracting 
the background (dashed line in (b)) fro m th e the difference 
function (data points in (b)), from Ref. 



locity space (V|| = Vj_= 0). This feature may correspond 
to the Coulomb ring associated with the proton emission 
from the heaviest fragment. In order to estimate quan- 
titatively the amount of evaporated protons from sec- 
ondary decay processes, one needs to estimate the con- 
tributions induced by other background effects that have 
to be properly subtracted. Few methods have been em- 
ployed to perform such backgrou nd estimation and sub- 
traction [lOl Eol ins lllll and they are all based 
on a scenario deduced from Boltzmann-Nordheim-Vlasov 
|l30j calculations. Multifragmentation is described as a 
two-step process. The first step is the cooling of the 
initial fused system through a sequence of light-particle- 
emission processes. The second step is the fragmenta- 
tion of the smaller remaining source where the remain- 
ing excitation energy is shared between a fixed number 
of primary fragments. These fragments then decay se- 
quentially while moving apart under the influence of the 
Coulomb force. An initial radial velocity can be added 
to the Coulomb motion in order to mimic a possible ex- 
pansion of the source. Based on such a scenario and on 
model predictions, two approaches have been suggested 
to estimate the background of the correlation function, 
differing in the way one compares calculations to exper- 
imental data. In one case, the shape of the background 
is deduced, normalized to the data and then subtracted 
[Tot. ,109. ,.110.. .lllj . in the other case, the input pa- 
rameters of the calculations are adjusted un til th e data 
points are reproduced (kind of back-tracing) |108| . After 
the background is subtracted, one can access the amount 
of evaporated particles from the primary fragments. In 
the following subsection we will describe how the tech- 
nique is applied to experimental data and what physical 
information has been extracted. 



B. Application to data and experimental results 

Figure 29 shows the experimental relative velocity cor- 
relation function (a) and the difference function (b) of 
the phosphorus - alpha pairs measured in the central col- 
lisions Xc + Sn at E/A = 32 MeV. To build the correla- 
tion function shown in this figure, Eq. (2) has been used 
by replacing the momenta pi and p2 with the velocities vi 
and V2 of a given fragment (here is a phosphorus) and the 
LCP (here an a particle). The uncorrelated events have 
been constructed by the event-mixing technique : For a 
given fragment in an event i having n LCP, we take ran- 
domly n LCP from n different events, then the relative 
velocity of this fake event is ca l cula t ed. The dotte d lines 
represent the background [lol [loi [lol [lll[m| . This 
background should contain any other contribution, which 
is not evaporated from parents of phosphorus. The rela- 
tive velocity distribution presented in Fig. 29c, obtained 
by subtracting the background from the difference func- 
tion, can be considered as the deduced contribution of a 
particles evaporated by parents of phosphorus. 

This procedure was applied to all fragment-LCP pairs 
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FIG. 30: Average secondary multiplicities per fragment of the 
evaporated hydrogen and helium isotopes as a function of the 
atomic number of the detect ed fr agments for Xe+Sn central 
collisions at E/A = 50 MeV fl07| . 



produced in the central collisions Xe + Sn at E/A = 32, 
39, 45 and 50 MeV. Thus, for each detected fragment, 
the average multiplicities of the light particles evaporated 
by the primary fragment was determined. These multi- 
plicities increase slightly with the size of the fragment. 
However they remain weak, not exceeding the value of 
1.5 at all incident energies, implying that the excitation 
energy of the primary fragments is particularly moder- 
ate. An example of the extracted multiplicity is given 
in Fig. 30 for central colhsion of Xe -I- Sn at E/A = 50 
MeV. The multiplicity of a given LCP does not change 
much with incident energy. From the spectra of the evap- 
orated light charged particles, the average kinetic energy 
has been extracted. 

Both observables (average multiplicity and kinetic en- 
ergy) have been used in order to deduce the average 
charge of the primary fragments, < Zpr >, their aver- 



age mass, < Apr > and their excitation energy, E*. The 
average charge of the primary fragment is calculated as 
the sum of the charge of the detected fragment, Z and 
the charge of all light charged particles correlated to that 
fragment, < Zlcp >■ This last quantity is defined by 
< Zlcp >= 'l2i^i' < Mi > , where Zi and < Mi > 
are the charge and the average multiplicity of the evap- 
orated particles, i = p, d,t,^ He, a. Since the isotopes 
of the detected fragments are not resolved and the neu- 
trons are not detected, two extreme assumptions were 
necessary in order to determine the primary mass of the 
fragment. The first assumption consists of assuming that 
the detected fragment is produced in the valley of stabil- 
ity. The second assumption considers two cases: i) the 
primary fragments are produced in the valley of stabil- 
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FIG. 31: Average excitation energy per nucleon of the pri- 
mary fragments as a function of their charge for central 
Xe-I-Sn collisions at E/A = 32, 39, 45 and 50 MeV. The hor- 
izontal lines represent a fit to the average excitation energies 
with a Z-independent value for each bombarding energy. Left 
panels: the primary fragments have the same N/Z as the com- 
bined system. Right panels: the fragments are produced in 
the valley of stability. The masses of detected fragments are 
assumed to follow the valley of stability except for the data 
points r epre sented by open circles where the EAL assumption 
is used lllf . 



ity; ii) the primary fragments keep the same ratio N/Z as 
the initial system. The obtained charge of the primary 
fragments varies between 1 and 5 units of charge in ad- 
dition to the charge of the detected fragment. The mass 
of the primary fragments depends also on the considered 
assumptions. The average multiplicity of the neutrons is 
deduced by the mass conservation, knowing the mass of 
the primary fragment, of the detected fragment and that 
of the secondary light particles. Multiplicities of about 
7 are reached, but depend strongly on the assumption 
made on the mass of the primary fragment. 

Finally with the help of all these variables, it was pos- 
sible to apply a calorimetry method to determine the ex- 
citation energy of the primary fragments E*^.. Figure 31 
shows the average excitation energy per nucleon E*j./A of 
the primary fragment as a function of its detected charge, 
and this for four incident energies and the two assump- 
tions on the mass of the primary fragment. The horizon- 
tal line in Fig. 31 represents the average value < e*^ > 
within the range of all studied primary fragments. Be- 
sides nuclei with low charges, all the experimental points 
lie on this straight horizontal line within error bars. In 
other words, whatever the incident energy and the mass 
assumption considered, the excitation energies per nu- 
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FIG. 32: Average excitation energy per nucleon of primary 
fragments produced in central collisions of the Xe+Sn system 
as a function of incident energy. The black and open symbols 
correspond to two hypotheses of the fragment masses. The 
vertical line s represent the standard deviation from the mean 
values im. 



cleon of the primary fragments are constant. Figure 32 
shows the evolution of the average value with the bom- 
bardment energy. The vertical bars represent the stan- 
dard deviations from the mean values. They are small 
and do not exceed 1 AMeV, consolidating the constancy 
of the value < e*^ >. In the case of the N/Z conservation 
assumption, the excitation energy per nucleon increases 
from 2.2 AMeV at E/A = 32 MeV up to a saturation 
value of 3 AMeV at E/A > 39 MeV. In the case of the 
other assumption, where the primary fragments are pro- 
duced on the valley of stability, the values saturate also 
but at lower energy. The constancy observed in Fig. 32 
of the excitation energy per nucleon for the various pri- 
mary fragments, can suggest that thermodynamical equi- 
librium was reached during the disintegration of the sys- 
tem. On the other hand, the saturation of < e*^ > can 
indicate that the fragments have reached their limiting 
excitat ion energ y per nucleon (or their limiting temper- 
ature) pllli3a. 

How to choose between the two assumptions of mass 
on the primary fragments? This choice was dictated by 
calculations with a statistical code GEMINI 133]. This 
code is very well suited at low excitation energies not ex- 
ceeding 3-4 AMeV. The procedure consists in using the 
deduced primary- fragment characteristics, Zpr, Apr and 
E*^ as input parameter, letting them decay and compar- 
ing the resulting multiplicities of the evaporated par ticles 
to the experimental values. It was shown in Ref. |l07j 
that the assumption of the N/Z conservation was reason- 
able. 

The average multiplicities of secondary light charged 
particles allowed an intrinsic characterization of the pri- 
mary fragments, but it can also be used to give valuable 
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FIG. 33; Atomic number of the heaviest fragment in the event 
as a function of center of mass parallel velocity, for Xe-|-Sn at 
E/A = 100 MeV reaction. 



information on multifragmentation events. Indeed, the 
ratio of the secondary evaporated LCP multiplicities to 
the total detected LCP multiplicities gives the fraction of 
therrnally produced particles. It has been shown in Refs. 
[lid llllj l that the maximum proportion of evaporated 
particles does not exceed, on average, 35% of the total 
number of produced light charged particles. For incident 
energies between E/A = 32 and 39 MeV, the proportion 
of evaporated particles increases, thus reflecting the in- 
crease of the excitation energy of the fragments observed 
in Fig. 32. Beyond E/A=39 MeV, this proportion de- 
creases to reach 23% at E/A=50 MeV while the excita- 
tion energy of the primary fragments remains constant. 
It should be noted that the extracted proportion of sec- 
ondarily evaporated particles constitutes a lower limit, 
because it does not contain the contributions which can 
come from the disintegration of the unstable nuclei like 
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FIG. 34: Evaporated proton multiplicity as function of the 
atomic number of the emitter: quasi-projectile or fragments 
formed in central collisions. The systems presented here are 
indicated in the figure. 



26 




FIG. 35; Comparison between SMM calculations (curves) and data (symbols) are presented. Left panel: Average excitation 
energy of the primary fragments as a function of their atomic number for central Xe+Sn collisions at E/A— 32-50 MeV. Right 
panel: Total charge co ntrib utions of secondary evaporated particles Zlcp- The dotted histogram represents the calculated 
freeze-out contribution 



the ^B e, ^L i etc. and the decay of the short-hved excited 
states |l34j . Similar results have been extracted also from 
the study of Kr-fNb reactions at E/A = 45 MeV [lo3| . 
It has been shown that i) the excitation energy of the 
primary fragments does not exceed 2.5 AMeV, ii) about 
80% of the detected LCPs do not originate from the sec- 
ondary statistical decay of the primary fragments. 



C. Comparison to models 

The described techniques have been applied to the 
quasi-projectile (QP) nuclei formed in Xe + Sn collisions 
at E/A = 100 MeV in order to extract their thermal 
characteristics |128| . Figure 33 gives an overview of the 
collision. It shows the atomic number of the heaviest 
fragments in the event as a function of its velocity par- 
allel to the beam direction. One can observe clearly two 
main components, one centered at projectile velocity and 
the other one centered at target velocity (the heavy quasi 
target fragments are not detected since their velocity do 
not exceed the energy threshold of the used detectors). 
Between the two components one can also observe an 
emission of fragments at mid-rapidity and characterized 
by a low charge and a high cross section. 

This preliminary result is surprising since the excita- 
tion energies of the spectator at E/A=100 MeV and par- 
ticipant fragments at E/A=50 MeV are almost the same. 
Is it an indication of thermal energy saturation? Or are 
we dealing with the same fragment production mecha- 
nism? The question is open. 

The correlation function method has been applied to 
extract the secondary evaporated light charged particles 
from each QP. Figure 34 shows results for two systems 



having different isospin: ^'^^Xe + ^^'^Sn, ^^'^Xc+^^'^Sn. 
The figure represents the extracted proton multiplicity 
as a function of the charge of the QP emitter. The data 
points obtained above for central collisions and same sys- 
tem but at a lower beam energy E/A=50 MeV are super- 
imposed on the same figure. All data points follow the 
same systematic: the evaporated proton multiplicity in- 
creases with the charge of the quasi-projectile or with the 
charge of the fragment (in the case of central collisions). 
The proton multiplicity values do not exceed 1.5 and are 
compatible with an excitation energy of 2-3 AMeV. 

The experimental estimate of the secondary decay 
component can be compared with the predictions of 
statistical multifragmentation models such as the SMM 
'l2l"| or the micro-canonical model of Raduta et al. 
ll22, 123, 124] which explicitly consider fragment exci- 
tation. The comparison of the extracted quantities with 
the calculations can indeed constitute a crucial test o f the 
basic assumptions of these models. In the MMMC |ll5l | 
approach, the excitation of primary fragments is implic- 
itly included in the neutron production at freeze-out, and 
a direct comparison is not possible. 

Statistical Multifragmentation model (SMM) calcula- 
tions have been performed in order to reproduce the 
thermal component. The version of SMM used in 
llldllllj gives access to the freeze-out configuration, i.e., 
to primary-fragment characteristics before secondary de- 
cay. Figure 35 shows the results of this calculation. The 
left panel in this figure represents the excitation energy 
in MeV of the primary fragments as a function of their 
atomic number (same data points as in Fig. 31 with the 
only difference being the use of MeV/nucleon as units of 
excitation energies). The right panel represents the total 
charge contributions of secondary evaporated particles, 
Zlcp- The data indicate that a maximum LCP evap- 
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FIG. 36: Proton and alpha energy spectra produced in cen- 
tral collisions of Xe+Sn at E/A=50 MeV. The symbols are 
the data and curves represent AMD calculations. The dots 
(continuous histogram) represent the total energy spectra, the 
triangles (dotted histogram) represent the secondary evapora- 
tion contribution and the open symbols (hatched histogram) 
represent the direct early emission. 



oration is obtained for E/A — 39 MeV and a decrease 
of thermal contribution takes place above this energy. 
Although the internal excitation energy of the primary 
fragments is well reproduced at each incident energy, the 
trend of the total charge of secondary LCPs, ZlcPi with 
the beam energy is not reproduced. This discrepancy 
can be understood if we consider the increasingly impor- 
tant effects of the collision dynamics as the beam energy 
increases. Direct emissions of LCP increase with increas- 
ing incident energy, while the proportion of the thermal 
contribution decreases. 

The obtained results have also been compared to dy- 
namical calculations. In order to understand the mecha- 
nism responsible for the saturation of the excitation en- 
ergy of the p rimary fragme nts observed in Fig. 33, AMD 
calculations |ll6l Ill7l IllSj for central collisions of Xe -I- 
Sn at E/A=50 MeV have been performed. A reasonably 
good agreement for the charge distribution, the charge of 
the heaviest fragment in the event and ave rage kinetic en - 
ergy of the fragments has been obtained jllOl Il28l Il35l| . 
The comparison of such calculations to the experimen- 
tal data is shown in Fig. 36. It represents the energy 
spectra of proton and alpha particles produced as free 
particles in AMD (direct emission) and secondary parti- 
cles evaporated from the calculated decay of the excited 
primary fragments. The two contributions are compared 
to experimentally deduced secondary particles with the 
method described above. AMD predicts reasonably well 
the total energy spectra of protons (78% of total emis- 
sion), but it fails to reproduce the thermal contribution 
(22%). 



XI. CONCLUSIONS 

The study of heavy-ion collisions at intermediate ener- 
gies has provided important information about the mech- 
anisms of multifragmentation and their links to a possi- 
ble nuclear liquid-gas phase transition. Our capability of 
fully understanding these phenomena depends strongly 
on how well one can identify and characterize the ther- 
mal and dynamical properties of fragment and particle 
emitting sources. In this chapter we have presented a 
review of the most significant results that different re- 
search groups have achieved in the last decades by using 
two-particle correlation techniques. 

Intensity interferometry techniques have demonstrated 
that we have space-time probes capable of measuring sizes 
as small as 10~^^m and time intervals as short as several 
10~^^sec. This can be viewed as an important advance 
in the field of heavy-ion physics. These space-time char- 
acterization techniques have provided valuable informa- 
tion about the size/volume, the density of decaying nu- 
clear systems, as well as quantitative estimates of emis- 
sion lifetimes and particle chronology. The extracted re- 
sults demonstrate that the emitting sources produced in 
heavy-ion collisions form a quite complex system: differ- 
ent particles are emitted at different times and by dif- 
ferent mechanisms. In this respect, it seems clear that a 
complete space-time characterization of nuclear reactions 
requires a study of multiple correlations extended to all 
particle species and including neutrons. Proton-neutron 
chronology is also indicated as an important candidate to 
explore the density dependence of the symmetry energy 
which represents one of the most challenging perspectives 
for the future. 

The recent introduction of imaging techniques has cer- 
tainly changed our interpretation of correlation observ- 
ables. The source sizes extracted with this technique can 
differ significantly from those extracted from other ap- 
proaches. This difference suggests that a detailed shape 
analysis of correlation observables is necessary in order 
to obtain a correct determination of nuclear densities 
and lifetimes. In this chapter we have also shown how 
techniques were found to quantify the strength of sec- 
ondary decays. Light particle-fragment correlations have 
been used to extract information about the excitation 
energy of primary unstable fragments. The possibility of 
constraining secondary decays represents an important 
perspective: several observables are indeed affected by 
the unknown contributions to the experimental spectra 
from different emission mechanisms dominating different 
stages of the reaction. 

Despite the large amount of information that the sci- 
entific community has been able to extract, certain as- 
pects have not been fully explored and require further 
research. For instance, from an experimental point of 
view, higher isotopic resolution could provide an im- 
portant opportunity to investigate isotopically resolved 
fragment-fragment correlation functions. This perspec- 
tive is very important in view of clarifying the indications 
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of evolutionary frcczc-out conditions in heavy-ion colli- 
sions and their possible links to the equation of state of 
asymmetric nuclear matter, relevant to astrophysical en- 
vironments like neutron stars and supernova explosions. 
Higher angular resolution and a large solid angle cover- 
age arc important requirements in order to increase the 
quality of measured correlation observables and to ex- 
plore their features while having a complete character- 
ization of the collision event (impact parameter, reac- 
tion plane, exact determination of velocity vectors, etc.). 
Higher resolution devices will be the key to the perspec- 
tive of extending imaging techniques to several particle 
species, improving complex particle correlation analyses 
that provide information about the space-time proper- 
ties at freeze-out and about the excitation energies of 
primary fragments. Other fields of nuclear physics re- 



search will certainly profit from the existence of detec- 
tor setups that are characterized by outstanding corre- 
lation capabilities. This is particularly the case for the 
nuclear structure groups working on spectroscopic prop- 
erties of exotic nuclear systems explored with the future 
radioactive-ion-beam accelerator facilities. 

The close interaction between the theoretical and ex- 
perimental communities will certainly contribute signifi- 
cantly in improving our capability of characterizing emit- 
ting sources. The use of more powerful experimental se- 
tups and the implementation of full quantum multi-body 
approaches to correlation functions promise to provide 
unambiguous information about the thermal and dynam- 
ical scenarios characterizing multifragmention phenom- 
ena and their links to the equation of state of nuclear 
matter. 
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